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ABSTRACT 

Estimates  for  the  derivatives  of  the  stresses  in  thin 
shells  and  plates  have  been  derived  by  Fritz  John  in  the 
case  where  there  are  no  body  forces,  no  surface  forces,  and 
the  medium  is  of  constant  thickness. 

The  procedure  used  has  been  modified  to  derive  similar 
estimates  in  the  case  of  a  continuum  of  variable  thickness 
acted  upon  by  body  and  surface  forces.   It  is  assumed  that 
the  strains  and  external  forces  on  the  undeformed  shell  are 
bounded  and  small,  and  that  the  curvature  and  variation  in 
thickness  of  the  shell  are  not  too  large.   As  a  result  solely 
of  these  assumptions  and  the  three-dimensional  equations  of 
non-linear  elasticity  rigorous  estimates  are  obtained  on  the 
stresses  and  their  derivatives  in  a  region  away  from  the 
lateral  edges  of  the  shell.   These  estimates  depend  only  on 
the  strain  energy  density  function  and  the  physical  dimensions 
of  the  shell. 

Starting  from  the  mathematical  three-dimensional  equations 
of  non-linear  elasticity,  one  can  derive,  as  a  consequence  of 
these  estimates,  a  set  of  two-dimensional  approximate  dif- 
ferential equations  for  shells  with  precise  estimates  on  the 
error.   These  new  equations  take  into  account  the  variable 
thickness  of  the  shell  and  loading  on  the  surface 
and  they  contain  as  a  consequence  those  of  v.  Karman  for 
plates  of  constant  thickness  with  normal  loads. 


Introduction 

We  consider  a  perfectly  elastic  homogeneous  and  iso- 
topic  shell  in  equilibrium  with  applied  surface,  edge,  and 
body  forces.   The  faces  of  the  shell  in  the  unstrained  state 
shall  be  of  variable  distance  from  the  middle  surface.   Under 
the  action  of  the  applied  forces,  the  shell  deforms  to  a  new, 
strained  position. 

We  apply  techniques  derived  from  the  theory  of  elliptic 
partial  differential  equations  to  obtain  a  priori  estimates 
for  stress  components  and  their  derivatives  which  .can  be  used 
to  formulate  equations  which  describe  the  deformed  state  of 
the  shell. 

These  estimates  have  been  established  by  Fritz  John  [3] 
for  shells  of  constant  thickness  with  no  applied  forces.   The 
techniques  used  are  extended  and  modified  to  establish  similar 
estimates  for  an  elastic  layer  of  variable  thickness  with 
applied  forces.   No  ad  hoc  assumptions  about  the  relative 
orders  of  magnitude  of  the  transverse  and  longitudinal  stresses 
and  their  derivatives  ( "Kirchof f  hypotheses")  are  employed. 

The  equations  formulated  reflect  the  variations  in 
thickness  of  the  shell  and  the  presence  of  body  forces.   These 
equations  are  similar  to  the  classical  equations  due  to 
V.  Karman  which  have  been  used  to  study  plates  of  variable 
thickness,  plates  on  elastic  foundations,  and  buckling  of 
plates.   In  the  case  where  the  shell  is  of  constant  thickness 
and  there  are  no  applied  forces,  the  formulated  equations  are 
similar  to  those  derived  by  Fritz  John. 


A  knowledge  of  the  geometry  of  the  shell  in  its 
unstrained  state  and  some  elastic  description  of  the  shell-- 
stress  strain  law--and  bounds  on  the  applied  forces  are  all 
that  are  needed  to  determine  these  estimates. 

If  the  thiclcness  of  the  shell  is  small  and  does  not 
vary  too  wildly,  the  curvature  of  the  middle  surface  is  not 
too  large,  and  the  strains  are  small  everywhere  in  the  shell, 
then  certain  estimates  on  the  stresses  and  their  derivatives 
are  valid. 

To  make  our  assumptions  rigorous,  we  define: 

h  =  undeformed  variable  distance  to  the  center  plane  of  the 

shell  from  the  surface; 
R  =  "generalized"  radius  of  curvature,  a  constant; 
D  =  distance  from  the  point   P  under  consideration  in  the 

shell  to  the  lateral  boundary  ("edge")  of  the  shell; 
H  =  maximum  of  h  over  the  region  of  consideration; 
H*  =  minimum  of  h  over  the  region  of  consideration; 
e  =  upper  bound  for  the  strain  everywhere  in  the  shell. 


The  quantity 


r2H    ZH    H  /-"I 
^  ^  V2H*R  " 


can  be  used  as  a  common  measure  for  the  "thinness,"  "shallow- 
ness" and  "variation  of  thickness"  of  the  shell,  and  of  the 
linearity  of  the  generalized  Hooke's  law  at  the  point   P  in 
the  undeformed  middle  surface  of  the  shell,  which  is  the  cen- 
ter of  a  disk  of  radius   D. 

If  9   is  small,  a  two-dimensional  approximate  description 


of  the  deformed  state  of  the  shell  becomes  feasible.   The 
quantity  X     defined  by 


y   =  -^  =   min(^j,  V2K*R,  J 


provides  a  lower  bound  for  the  "wave  length"  of  the  stresses 
and  strains  in  the  longitudinal  directions. 

In  general  it  will  be  shown  that  the  n-order  longitu- 
dinal derivatives  of  the  strain  components  are  bounded  by 
K  e/X*^  when  9   <  9_,   and  the  constants   9^,  K,  ,  K^,  ••• 
depend  on  the  elastic  properties  of  the  material  alone.   In 
addition,  conditions  will  be  imposed  on  the  order  of  magnitude 
of  the  applied  forces  in  order  that  the  stresses  internal  to 
the  shell  be  of  similar  order.   If,  as  assumed,  the  maximiim 
stress  is  of  order  E   (where   E  is  Young's  modulus),  it 
can  be  shown  that  the  transverse  shear  stress  are  of  order 
Eey  9t  each  point  in  the  interior  of  the  shell,  and  that  the 
longitudinal  stress  varies  linearly  within  an  order  of  error 

of  magnitude   Es— ^  along  a  normal  fiber  of  the  shell.   Here 

X 
the  expression  "order  of  error"  indicates  that  any  constants 

tacitly  involved  in  the  estimate  can  at  most  depend  on  the 
material  (stress--strain  description  chosen). 

Once  these  estimates  have  been  obtained,  approximate 
two-dimensional  equations  ("interior  shell  equations") 
involving  surface  and  body  forces  can  be  derived  from  the  exact 
equations  of  a  three-dimensional  elastic  media.   This  is  accom- 
plished by  grouping  terms  in  the  equations  according  to  powers 


By  wavelength  of  a  function  f ,   we  mean  a  quantity  propor- 
tional to  max|f |/|f • I . 


of  h/x  appearing  in  the  concrete  estimates,  and  hence  the 
errors  of  approximation  will  be  known. 

The  technique  used  is  that  of  deriving  bounds  on  the 
Hilbert  norms  of  the  derivatives  of  the  stresses  from  energy- 
identities.   These  bounds  can  be  converted  into  pointwise 
estimates  by  means  of  a  Sobolev  type  lemma.   Since  the  equa- 
tions involved  are  non-linear,  it  is  clear  that  estimates  for 
the  norms  must  involve  this  kind  of  technique. 

In  proposing  a  system  of  shell  equations  we  choose  as 
the  unknown  functions  the  three-longitudinal  components  of 
stress  in  the  middle  surface  and  the  three  coefficients  of 
the  second  fundamental  form  of  the  deformed  middle  surface. 
The  three-dimensional  theory  of  elasticity  furnishes  six  non- 
trivial  approximate  equations  for  these  six  quantities. 

If  in  addition  the  shell  is  very  shallow,  then  these  six 
approximate  equations  can  be  further  reduced  to  two  equations 
for  two  unknown  functions.   One  of  these  is  analogous  to  the 
Airy  stress  function,  and  the  other  is  related  to  an  approxi- 
mate integral  of  the  second  fundamental  form  of  the  deformed 
middle  surface  of  the  shell.   In  the  case  where  the  shell  is 
of  constant  thickness  end  has  no  initial  curvature,  this 
system  of  equations  agrees  with  that  ascribed  to  v.  Karman  [4] 
and  F3ppl.   In  all  the  systems  of  equations  mentioned  above, 
terms  are  included  which  arise  from  the  variable  thickness  of 
the  shell,  and  the  body  and  surface  forces.   In  particular,  a 
term  due  to  normal  equal  and  opposite  loading  of  the  shell 
appears. 


Chapter  1  -  Summary  of  the  Equations  of  a  Homogeneous 
Isotropic  Perfectly  Elastic  Solid  in  Equilibrium, 
with  Body  and  Surface  Forces 


1.1  Coordinate  System 

We  use  Lagrange  coordinates  U  =  (U  ,U  ,11"^)   to  dis- 
tinguish the  individual  particles  of  the  solid.   The  solid  is 
considered  in  two  states — "unstrained"  and  "strained."  Let 
the  Cartesian  coordinates  of  the  same  particle  be  X  =  X  (u) 
for  the  unstrained  state,  and  x  =  x"^(u)   for  the  strained 
state.   By  f  .   we  denote  the  partial  derivative  of  the 
function  f  with  respect  to  the  variable  U  .   Line  elements 
in  the  unstrained  and  strained  states  respectively  are 

dS^  =  dX^dX^  =  Gj^j^dU^dU^  and 
ds^  =  dx^dx^  =  g^j^dU^dU^,   where 

'  ^ik  =  x',i^^k  ^^^  gik  =  ^',i^',k 

are  the  unstrained  and  strained  metric  tensors.   Although 

tensor  operators  can  be  performed  in  either  metric,  all 

shifting  of  indices  will  be  performed  with  the  unstrained 

metric  tensor  G.,  ,   which  in  principle  is  known.   In  addition 

ik 
2 
the  dS   metric  can  be  used  to  express  covariant  differenti- 

ation  with  respect  to  the  variables   U^   in  the  ds   metric. 

2 
The  Christoffel  symbol  for  the  dS   metric  is 


-See  John  [3],  Truesdell  [7]. 


k 


2 
and   for  the     ds        metric 

_1  1— ir/  \    _     i  i 

^jk  =   2^      ^^rj,k  "^  ^r\^,j    ■   Sjk,r^    =  ^Jk  +  ^jk 

where  c  .,   is  the  tensor: 

^^^  ^jk  "  2^  ^^rj;k  "^  ^rk;j  "  ^kj;r^' 

";k"   denotes  covariant  differentiation  with  respect  to  U 
in  the   dS   metric,  and   g"""^  and   g  ^     are  defined  by 

/o„^  — ir      c-i    i^   -im^rs 

(2a)  g     &^^=   \,      g        =  G  G  g^3. 

The  conditions  of  compatibility  reflect  the  fact  that 

2         2 
ds   and  dS   are  Euclidian  line  elements.   Expressing  the 

2 
vanishing  of  the  RLemannian  curvature  tensor  for  the  ds 

2 

metric  in  terms  of  the  dS   metric,  gives  six  independent 

conditions : 

.      1  .abba. 

O)      0  =  2^Shk;ij'*'Sij;hk"2hj;ik"^ik;hJ^'^^ab^'^iJ^hk"^hJ^ik^' 

The  principal  strains  are  the  eigenvalues  of  the  strain 
matrix  defined  by 


1.2  Strain  Energy  Density  Function  W, 
Stresses,  and  Equations 

For  a  homogeneous  isotropic  perfectly  elastic  material, 
the  strain  energy  density  per  unit  unstrained  volume  is  a 




Numbers  such  as  J>.2   will  refer  to  Chapter  3>  item  2 


symmetric  function  of  the  principal  strains.   We  require  that 
W  agree  with  the  classical  strain  energy  density  function  for 

small  strains. 

rnk 
The  contravariant  components  of  the  stress  tensor  t 

in  the  U   coordinate  system  are  defined  in  such  a  way  that 

the  force  acting  on  an  element  of  area  ds  with  direction 

cosines   5   in  the  strained  state  has  Cartesian  x-components : 

(5)  *"'^^\m'''',k?r'^^'   ^°^  ^  =  ^'^'^• 

The  stress-strain  relations  are  given  by 

with   (g/G)  =  (det  g^j/det  G^^.)  =  l+2s^+2  (s^-S2)+4(2s^-3s^S2+'^sp 

_i      _ik      _ijk 
^1  "  ^i*  ^2  ~  ^k^i'   ^3  ~  ^j^k^i- 

The  equations  of  equilibrium  take  the  form 

where   ";;k"  denotes  the  covariant  derivative  with  respect  to 
u   formed  in  the  ds   metric,  and   eP  is  the  u"^- component 
of  the  body  force  per  unit  deformed  volume  (normalized  to 
agree  with  the  unit  of  force  chosen  for  the  stresses). 


1.3  Differential  Equations  of  Elasticity 
in  Pseudo-classical  Form 

The  tensor  or  invariant  form  of  the  differential  equa- 
tions are  inconvenient  to  use  in  estimating  the  derivatives  of 


the  stresses.   Instead  we  regard  the  equations  as  between 

ik 
vectors.   An  object  such  as   t  =  (t   ),   with  any  number  of 

subscripts  or  superscripts,  can  be  regarded  as  a  vector — a 

point  in  a  normed  linear  space  of  dimension  three  times  the 

number  of  indices    The  norm  of  a  vector  shall  be  the  square 


root  of  the  sum  of  the  squares  of  its  components.   For  example, 

ik,    _._ ^  ..   _   ,  ,    /..ik  ik 

rs' 


if  w  =  (w  „),   then  the  norm  of  w  is   Iwl  =/w^^w^ 
rs '  '     '     \/      vs    2 


ik 
Vectors  can  be  multiplied  by  scalars:   Xw  =  ( Xw   ) .   Vectors 

of  the  same  type  can  be  added:   if  u  =  (u„_)   and  v  =  (v„_) 

rs  rs 

ik   ik 
then  u  +  V  =  (u  +v   ) .   Any  two  vectors  can  be  multiplied : 
^  rs   rs  ^ 

/  ik.  ,  ab,   .         ,  ik  ab. 

if  w  =  (w   )   and   v  =  (v   )   then  wv  =  (w  v   ), 
^  rs  c  ^  rs  c 

2        /  ab  de .                    2    /  ik  An. 
V  =  (v  V-  )   and   w  =  (w  w   ).   The  gradient  of  a  vector 
^  c   f  rs  np'       ^ 

w  with  respect  to  the  Lagrange  coordinates   U   is  again  a 

ik  ik 

vector,  denoted  by  w' :   if  w  =  w   ,   then  w'  =  w   . .   We 

rs  X s  J  w 

also  have  for  any  two  vectors  u  and  v:   |uv(  =  |u|»|v|. 

In  many  cases,  only  the  general  form  of  an  expression 
will  be  of  interest.   We  use  the  following  notation  to  dis- 
play this.   If  u,  V,  w,  p,  q  are  vectors  then  F(p,q) [u+v+w] 
will  denote  a  vector  valued  function  each  of  whose  components 
is  a  linear  form  in  u,  v,   and  w,   with  coefficients  that 
are  functions  of  p  and   q,   defined  and  differentiable  as 
often  as  needed  for  sufficiently  small  p  and  q.   F  may, 
and  in  general  will,  stand  for  a  different  expression  in  each 
equation. 

The  assumptions  on  W  indicate  that  its  form  is  for 
-1<  V  <  I 

8 


Adjusting  the  units  of  force  normalizes  Young's  modulus   E, 
leaving  Poisson's  ratio  as  the  only  classical  elastic  constant. 

Introducing  t  =  (t,  ),   and   e  =  (Sj^),   we  have  by 
(1.8),  (1.6), 

(9)  <=(l-2v)(l.v)^;-I^i-^(-)-^' 

Solving  for  e   in  terms  of  t  we  obtain  a  stress-strain 
relationship  of  the  form: 

(10)  e^  =  (l+v)t^  -  vtj5^  +  F(t)t^. 

,  i. 
We  introduce  the  vector  t^  =  (^  )   by  the  highly  non-invariant 

equation 

(11)  G^k  =  \   +  ^k- 
Then, 

(12)  G^^  =  \  -  ^k  "^  F(Tl,t)Tit,   r^j  =  F(ti)ii', 

(13)  t^^  =  *k  "^  F(n,t)Tlt  =  t^^    +    F(Tl,t)T]t, 

and  from  (1.4)  and  (l.lO)  we  have 

(14)  gik  =  ^ik"^^ik  =  G^k'^2(l+v)t^j^-2vt^.j5j^+F(Tl,t)[t%t], 
and  from  (1.2),  (1.12),  (l.l4) 

(15)  c^^  =  P(n,t)[t'+Ti't]. 

The  equations  of  equilibrium  (1.?)  become  using  (1.15),  (1.12) 

(16)  t^'"   +  B^  =  F(Ti,t)[tt'+Ti't+Tit»]  =  P^. 

,  m 


The  conditions  of  compatibility  (1.7)  yield,  upon  substitu- 
tion of  (1.11),  setting  i  =  j   and  use  of  (1.13),  (l.l6)  : 

=  F(Tl*t)[Tlt"+(t')^  +  Tl"t  +  (Tl')^t  +  Tl't'  +  tt"], 

which  resembles  the  Beltrami-Michell  compatibility  equations. 
In  each  of  these  equations,  the  expression  F(Ti,t)   stands  for 
a  different  function  uniquely  determined  and  regular  in  a 
neighborhood  of  the  Origin  which  is  dependent  solely  on  the 
choice  of  the  strain  energy  density  function  W.   We  remark 
that  by  (1.13),  it  makes  no  difference  if  t  =  (t^^)   or 
t  =  (t..  )   in  the  expressions  for  P   and  Q 

1 K 


lA      "(y"  and   "g"   Notation 

The  symbol  "i^"     will  be  used  in  two  ways.   First  in 
the  conventional  sense,  and  second  to  indicate  dependence 
only  on  W. 

The  equation 

A  =  [^(B),   B  >  0 

means  that  there  exists  a  constant  K  dependent  at  most  on 
the  choice  of  W  such  that   |A|  <  KB. 

The  symbol   "o"  will  be  used  only  in  combination  with 
"  ff  "  The  equation 

A  =  &(B)  +  o(C),   B  >  0  and   C  >  0, 


10 


means  that  there  exists  a  function  K(k)   defined  for  all 
positive  k  and  dependent  only  on  W  such  that 

|A|  <  K(k)3  +  kC   for  all  k  >  0. 

For  example, 

xy  =  (j  (|x|  )  +  a(|y|  )  where  a  possible  choice  for 
K(k)  =  -jTr-  independent  of  the  choice  of  W. 

Rules  for  the  use  of  the  symbols   "Pp"  ^"^^      "o"   easily 
proved  are: 

A  =  e^(B)  +  a(|A|)  -.  A  =  0^(8), 

A  =  (5(B)  +  a(C+|A|)  -^  A  =  &  (  B)  +  a(C), 

A  =  (5(B)  +  a(C)  -  A  =  &(&fC), 

A  =  ^(B)  +  a(C),   A'  =  (5'(B')  +  ct(C')  -  A+A'  =  8^(&hB')  +  a(C+C') 

A  =  0^(B)  +  a(C),   C  =  (^(D)  +  a(E)  -  A  =  0^(b)  +  a(l>4-E), 

A  =  (5(b)  +  a(C),   B  =  5(P)  +  o(E)  -  a  =  ^(d)  +  a(C+E). 


11 


Chapter  2  -  Description  of  the  Shell 

Greek  letters  a,  p,  7,  •••  will  be  used  for  indices 
ranging  over  1  and  2,  while  Roman  letters  i,  J,  k,  ••• 
will  be  used  for  indices  ranging  over  1,  2,   and  3- 

In  the  unstrained  state  the  boundary  of  the  shell  shall 

consist  of  two  faces  ^v, ,  ,   and  2    and  a  lateral  surface 

h+        h- 

("edge")   L.   We  assume  that  the  functions  h+  and  h-  are 
well  enough  behaved  so  that  it  is  possible  to  define  a  unique 
midsurface  2   of  the  shell  that  is  equidistant  from  the 
faces  21  _  and  ^v,,  •   We  consider  the  edge  L  to  be  composed 
of  normals  to  2  . 

We  choose  the  Lagrange  coordinates  U   by  the  condi- 
tion that  all  particles  on  the  same  normal  to  2   in  the 
unstrained  state  have  the  same  U^'  coordinates,  and  distance 
U   from  the  midsurface  2 

The  coefficients   E  „   and  L  „   of  the  first  and  second 
fundamental  forms  of  the  undeformed  midsurface  2:   will 
determine  2   to  within  a  rigid  motion.   The  metric  tensor 
of  the  unstrained  space  line  element  can,  as  a  result  of  our 
special  assumptions  on  our  "chosen"  Lagrange  coordinate 
system,  be  written  as 

Selecting  any  point   P_  of  ^^     as  the  origin  of  a 
Cartesian  coordinate  system  X   with  X   axis  normal  to  Z 
at  Pq,   we  assume  that  we  can  represent  the  surface  ^       by 
the  equation  X'^  =  f(X^,X  )   in  a  neighborhood  of  Pq.   By  a 
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rigid  motion  of  the  undeformed  Ca^-tesian  coordinate  system, 
we  can  arrange  that   f(0,0)  =  f  ^(0,0)  =  0.   We  make  the 
further  assumption  that  actually  the  entire  surface   2^  can 
be  represented  by  X^  =  f{X^,X^)     with  X°  varying  over  some 
domain  D* . 

We  define  the  thickness  of  the  shell — the  distance 

between  2,    and  2.    measured  perpendicular  to  the  mid- 

1   2 
surface  21  --as   2h(X  ,X  ).   Half  the  maximum  and  minimum 

thickness  will  be  defined  as 

H*  =   inf   h{X^,X^) ,      H=   sup   h(X^,X2). 
x"  in  D*  X^  in  D* 

We  associate  a  "characteristic  radius"   R  with  the 

surfaces  2   2   Z_   as  follows:   for  any  choice  of  the 

origin  P_  of  ^^,      X°^  in  D*  the  k  -order  derivatives 

of  f  with  respect  to  X*^  satisfy  everywhere  in  D*. 

(2a)        If^^^l  <  R-^(hR)^-^/^   k  =  2,3,^---. 

At  some  point  in  D*  such  that  X*^X^  <  2H*R 


(2b) 


r(f-h) 


<  (h/R)^/^ 


ax" 

Conditions  (2a)  and  (2b)  amount  to  the  fact  that  R  is  an 
upper  bound  for  the  radius  of  curvature  of  the  midsurface  and 
faces  of  the  shell,  and  that  at  some  point  in  the  shell  2 
and  2   are  approximately  parallel. 

Choosing  11°''  which  coincide  locally  with  X^'  we  have 
for  the  first  and  second  fundamental  forms  of  the  undefonned 

This  is  more  than  is  needed.   See  John  [3],  p.  244,  footnote  13. 
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shell  in  terms  of  f 


(u\u^) 


E„a  =  5°  +  f  ^f  ^, 

/■2\  op    p     ,a.    ,p 

L  ^  =  (1+f   f   )"-'-^^f  «. 

On  the  faces  of  the  shell  we  have  in  terms  of  %. — 
the  direction  cosines  of  the  normal  to  deformed  surface 

0  =  ?,dx^  =  kAx^     dU^  +  x^  ,h  du"). 
i       i^  ,CL  ,3  ,a   ^ 

Denoting  the  surface  forces  per  unit  deformed  surface 
area  acting  in  the  strained  state  in  the  x  coordinate  system 
by   S    on  2,   and  by  S    on  ^_,       we  have  from  (1.5) 
and  (2.3) 

S^-'  =    (t'"^-^  -  h  t"^°-^)x^-'  „(x^  ^U^, 
X    ^        ,a    '    ,m'   ,3^r' 

s^"  =  (t'^^-  -f  h   t"'«-)xi-  „(x^  ^^r, 

X    ^        ,a    '    ,m^   ,3  r' 
for  i  =  1,2,3. 

7s\]  1  + 

Multiplying  by  ^^     and  summing,  we  have  for  S"^    , 

S"^',   the  surface  forces  per  unit  deformed  surface  area 
acting  in  the  strained  state,  but  specified  in  the  Lagrange 
coordinate  system 

sJ-^=  (tJ>  -\a^^''"')(^',3^r)^ 
SJ-  =  (tJ5-  -^  h^^tJ°'-)(x^^^er)-. 
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Chapter  5  -  Parameters  Used  in  Making  A  Priori 
Estimates  on  the  Derivatives  of  the  Stresses 


Let   P-  be  any  point  on  the  mldsurface  ^   ,      and   D 
the  closest  distance  of  P^  from  the  lateral  surface  of  the 
shell  L.   Let   R  denote  the  characteristic  radius  associated 
with  the  middle  surface  2  .   We  have  by  our  assumptions  on 
h  that  H*  <  |h(U''",U^)  |  <  H,   where  the  shell  has  thickness   2h. 

Consider  strains,  shell  thicknesses,  surface  loads  S''' 

3  i-  3 

on  the  face  U  =  h  and  S    on  the  face  U  =  -h  which  are 

specified  in  our  Lagrange  coordinate  system  and  body  forces 

B  ,   also  specified  in  our  Lagrange  coordinate  system.   In 

the  region  M  =  { ( U""" , U^ , U^ )  | U^'u"  <  D^,  |U^|  <  hlU"*-,!;^)}   we 

require  that  these  quantities  satisfy  the  following  conditions 

given  in  terms  of  a  constant  e    ,      to  be  specified  in  the 

sequel: 

H^^i„i±        I  ,  „2„-m   ,  H\^,„i  ,  ,  2,,-m 


1    m  1    m-1 

where  i,k  =  1,2,3  and  m  =  0,1,.-.,N;   N  will  be  picked  as 
high  as  needed  in  the  sequel. 
We  define 

(2)  .  6  =  maxM,      -^,  ^yT). 

If  H,  H*,  R,  D  are  fixed,  then  0  is  a  non-decreasing 
function  of   e.   Let   e  be  determined  as  the  smallest  positive 
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number  such  that' 


,a± 


(3)   spec  ie^)    <   e,    |s"^^^ 

Is 


m 


,a,  •  •  •  a   —     0' 

1    m 

'  1    m 

'(e/9  )'""*' ■'■H""^e   if  the  number  of  U^   deriv- 
atives is  even  and  1/3,      or  the  number  of 
U   derivatives  is  odd  and  1=3, 


B 


»k, 


m-1 


<< 


(e/9Q)^'^^H~'^(h/H)e   if  there  is  no  U^ 


derivative  and   i  = 


_  1 


(9/eQ)'^H"  e(H/h)   otherwise, 


n  the  region  M,   and   m  =  0,1,  •  •  ■  ,li. 


H 


From  O.l   -   3-3)  it  follows  that  -^yi  <  9^     and  hence 

he 
0"  ^-""-"B  "  -  He, 


that   e  <  e..   Setting   n  =  rj-^-,   we  have  Q  <  -^  <   1 
We  define  the  wave  length  X,      and  ratio  Cl     as 


(M 


X  =  ^  =  e  min(D,  V/21FR,  ^,   "  =  x- 


Assumptions  (3.1)  imply  that 

(5)       spectrum(e,)  <e,   s""^     .     =  U  (^t)  » 

*^  '  1    m      X 

c3+  o3- 


'    1  m      X 


We  follow  Sensenig  [61 


16 


''       ^~¥^      ^^     i  =  3     and   the  number  of     U^ 

A 

differentiations  is  odd,  or  if  i  j^  3  and 


^  k    k    =  < 


the  number  of  U   differentiations  is  even 

(— )   if  there  is  no  differentiation  with 
^  m' 

respect  to  U   and  1  =  J> , 

( — ^ — )   otherwise, 
X   h 


for  m  =  0,1, • • • ,N. 

Restricting  9^     by  ^n  '^  2"  ^^  have 


(6) 


2H 


<  X  <  y^>      8hA  <  r,   X  <  D/2 


H* 


H*  1  .  n 
inx  -  x^ 


^  <  -^  <  -TT^T-  <  T»   ^"^d   £ 


.  (2H*)^  .  (2H)^ 


Since  f  and  its  first  derivatives  vanish  at   Pq,   we  have 

in.  a  region  MrU^'u^  <  X^  by  (2.3a)   |f  ^|  <v^x/R.   By 

(2.3c)    and   the   inequality      |  |h  „|-|f  „1|    <    ih  ^-f  ^1    <   vTTTR 

we  have      |h      |    <    v^i/R  +    \/2\/'R     at   some  point   in  the   region 
,  ot 

U°^°  <  2H*R  <  x^.   Since  by  (2.3b)   |h  „^|  <  R""^,   we  have 


,a3 


(7) 


|h  „|    =   QWW^  +   X/R)   =    0(^9). 


By   (2.3)   we  have  employing   (3.6),    (3-7) 


I  CX|  2/2 

l^ag    -    %^    <  2xVR   , 

IL^I    <  1/R. 


In  addition  we  have 


Eae,Y  =  Cr(x/R^),    E 


'aa>75 


=  (5^(1/R/ER),    E  „      _      =  (5(l/RhR),..., 


a6,75n 


^aa,T   =   0(1/R/HR),    L^^^^   =  (STd/RhR),.... 
We  have   from   (2.1)    and    (2.3) 
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(8)        n  =  0(h/R),      ^'    =  (5(1/R),      n"   =    &(1/F/H1),..., 
n*  =  0(h/}v1iR) ,     n"  =  6r(h/RhR) , . . . . 


If      ^Q  <  ^     we  have  by    (2.1) 


G  =    (det   E      )(l-U^L°'+(U^)^det   L   )^    >  0, 

ik 
hence      G  can  be   estimated  by 

G^^   -    5^  =&(|ti|)    =  &(h/R). 

For  the  Christoffel  vectors  r  =  (P..  ),   we  have 

(9)  r  =  gr(i/R),  r'  =  G^(i/r/Hr),  r"  =  @^(i/RhR),.... 

Our  assumption  of  the  smallness  of  the  principal  strains 
furnishes  the  estimate   e.  =  Kji&r^)      for  the  individual  elements 
e  ,   provided   9   is  sufficiently  small.   We  have  that 

K  0 

El,  =  G  e  ,  ,   G    positive  definite,   e  ,  =  e,  .   It  is  suf- 
k       mk        ^  '   mk    km 

i 
ficient  to  show  that   e    is  similar  to  a  symmetric  matrix. 

The  similarity  transformation  is  given  by  the  symmetric  posi- 
tive matrix  whose  square  equals   G   .   For   9   sufficiently 
small,   G    is  close  to   6.,   and  such  a  transformation  can 
be  found . 

The  strain  energy  density  function  W(s  ,s^,s^)   which 
we  consider  to  be  as  regular  as  needed  for  sufficiently 
small  |s  I --the  absolute  value  of  the  sum  of  the  i-powers  of 
the  principal  strains--determines  the  stress  strain  relations. 
It  follows  that 


\-^{^^4^ 
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Hence,  given  the  function  W,   bounds  t_^,  t\^     such  that 
|t|  <  t^  and   In  I  <  -n^  can  be  obtained  and  (I.9) ,  .^ .  .  ,  (l .  1?) 
are  valid  provided  8^     is  sufficiently  small. 
We  may  also  conclude  that  the  region 

M  =  t(U^,U^,U^)|U°U°'  <  X^,  lU^I  <  h] 

lies  entirely  within  the  shell.   The  distance  between  any 
point  on  2-  and   P_  is  given  by 

If  we  restrict  ^n  <  ir    we  have 

A/iA;°'+f2  <  -^ 

Hence  for  any  point  in  M,   the  distance  to  P„  is  at  most 


W  "^  T/^  ^  °'   provided   0^  <  ^. 


°  *h< 


In  what'  follows,  we  shall  use   e  to  denote  e_,   since 
no  confusion  will  result. 
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Chapter  k  -   Estimates  for  the  L„   Norms  of  the 
Second  Order  Derivatives  of  the  Stresses 


4.1  Coordinate  Transformation  to  a  Shell 
of  "Constant"  Thickness 

We  consider  the  transformation  of  the  Lagrange  coordi- 
nates U  =  (U"^,U^,U^)   to  the  new  coordinates   U  =  (ir",U  ,\r) 
given  by 

(1)  f  =  u^ 

h(U-',U^) 

Under  this  transformation,  the  center  plane  of  the  shell  2^ 

remains  fixed  and  the  faces  U  =  ih(U  ,u  )   are  transformed 

to  faces  2;    a  constant  distance  H  from  2  . 
±n  U 

ik 
For  the  moment,  we  choose  to  regard  the  t    as  tensors, 

''ik  i 

and  calculate  their  new  components   t    in  the  U   coordi- 
nate system.   From  the  transformation  formulas 

-ik   ^it    ^li,  im    .   ^An   au  3lf;ik 

t   =  n  — — t    and   t   =  — ^ rrt   ,   we  have 

air  au^  3U^  h\r 


(2)   t°^  =  t°^. 


.3a   h(u\u^)   3,      U^^ 


We  define 
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and  note  that  -^  =  -—-.      Consequently,  (4.2)  takes  the  form 


SU°    3U^ 


t°^  =  t^, 

(3)        t^°  =  t^^  +  (^6^)i^^ 

where   • ,a  = 


t--  =  t^^  +  (2h  «u)tP^  +  (h  u")tP^, 
—        ,p_—         >ct^ 


Using  the  chain  rule  for  differentiation  we  have 

Denoting  differentiation  with  respect  to  u"  by  the  symbol 
",a"  and  the  operator  ^  -^     by   ",3"  we  write  (4.4)  as 

(5)  •,k=  •,k  -  •,3^k(h,e^)'    -,3  =  M' 

For  the   Laplace   operator 

,    ..    =    .  +.,,=    .  +    .      -,,      where 

,ii  ,aa  ,33  ,aa  ,33 

•p 

=    .  -    .    ,    (2h     u)    +    .    ,,(h     h     u    )    -    .    -,(h u). 

,aa  ,00.  ,3o^      ,a-'  ,33      ,a   ,a-  '  ,3      >ga- 

We  remark  that  operators   ",3"  and   "ja"  are  not  commutative. 

Transforming  equations  (l.l6),  where  t    transforms  as 
a  contravariant  tensor  and  the  derivatives  transform  according 
to  the  chain  rule,  we  have  in  our  notation: 

(6)   t°'^^^  =  t°^^j^  =  -B^  +  P°, 

t^^  ^  =  t^^  ^  -  2h   (i)t^°  +  V  +  (-B^+P^)  =  -B^  +  P^, 
,K   —  ,K     ,an 

(7) 


V 


=  \a\i^  *   \ae.Ht*  -  h,„u(-B°.P°)  -  (-b'.P). 
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Using  (1.17),  the  "Beltrami -Michell"  equations,  and  the  fact  that 

t^^  -  (h  u)tP^  =  t^^  +  (h  u)t°^ 

}  P~~       "~       J  P. 

we  have 

(8)   t^^  ..  +  -—-{t^^   ,^  -  (h  u)t^^   )  =  J  = 

4^(^,^-)i^\ii-4''  -  ^k(^a^)tl^^ii^(^p.Hi^'),ii^J2^  ^ 

^k(^aii)^S3  -  ^3  -  ^^p^^(%3  -  ^05^^^' 
^^^^^  B^k  =  i^(Bj,j^^hk  ^  ^^k  ^  B^,h'  ^-^^  h  ^°^^^i^ 

at  most  first  order  derivatives  of  the  stresses  and  are  given 
by 


j^  =  (^5ii^i^'^i  ^  (^5^^ii^ 


53 


Derivatives  of  the  stresses  fall  into  two  classes  which 
are  dependent  on  the  number  of  3's  appearing  among  the  super- 
indices  and  subindices.   We  denote  the  N-order  derivative  of 
the  stress  vector  t   „    „   by  t^^^^^'   where   u  =  num- 
ber  of  3's  among  superindices,   s  =  number  of  3's  among  the 
differentiations,  and  define 
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1^(0)  (s)  ^  ^T^a3  ^ap  jl/2^ 

jN     ^   ,5 .  .  . ^ Og_^2^ .  .  .  Oj^   ,5 •  •  •  56  ,  .  .  .  o„ 

/Q>i      .(i)(s)  _  /-.5e  t^p  a/2 

t(2)(s)  ^  (^^33  ^  ^33  )l/2^ 

vrhere   u  =  0,1  or  2   and   0  <  s  <  N,   and  we  sum  over  all 
free  indices.   For  a^ny  vector  w  =  {\r  **'  )   we  can 

define,  (where  [i  =   number  of  3's  appearing  in  superscript) 

„(u)(s)  _  /V-33.  .  .3ap.  ..6 

^3...3a3...6  n1/2 

,3. .  -35^^^.  .  .5j^ 


and 


(u)   _  V,.,(u)(s)    e  _    V     ,Ju)(s) 


w^^'  ,,  =   )  w^"'\/,   w'..  =    )      w    ^ 


=  Yw(^^(^^   w^   = 

,N    L^        ,N  '   '',N 

s=0  U-s>0 

or  u-s  even 

Using  the  imles  of  transformation  given  by  (^'3)  we  have 

ik 
for  the  second  order  contravariant  tensor   t   : 

(10)  t(-)  =  Q(  I   aV^-^'^)  =  0^(t(^^  +  I   A'V^--^')). 

j=0  J=l 

We  note  that  by  (3-Da) 

(11)  A  =  (h^^h^gu^)^/-  =  &(fie),  h^g/h  =  tr(x"^), 

h^^u  =  &(Qe/x). 

We  define  for  a  multi  index  n,...!!-.  =  n  a  function 

1    N 

p. (N),   i  =  1,...,N!   the  set  of  all  possible  permutations 

N 
of  the  multi  index  n.   Obviously  for  any  function  w  e  C 
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«!n'=0(w!j|,„)  =  &(«!^>) 


for  s  =  0, . . . ,N,   and 


i  =  1,...,N!,   since  differentiation  in  the   U   coordinate 

system  is  commutative. 

a        "^ 
Differentiation  with  respect  to  U   and  \y      is  not 

2 
commutative.   We  have,  using  the  chain  rule  for  w  e  C  : 

w<|)=&(w(|)), 

w(o)=G-(w||).Gw;i');?'.<w|£))(^).<.wa))w 

Noting,  however,  that 


"!?)!i'=&a«!iO!?'-!i'«!iO- 


we  have 


„(o).g(w(|).(.wW)(°)..w(|).<..wU))(^)). 

For  w^_'   we  have  by  a  similar  calculation 
These  results  can  be  summarized  by 

~   i=i    ~        ^    ^1    ~ 

for   i  =  1,2   and   s  =  0,1,2. 

Combining  the  results  of  (^.5)  and  (4.10)  we  have 

(12)   t<M)(=^)  ^^(tMM    ,  Y  I     {A^V^-J')(^=') 

'  i    £=0  J  =  5° 
where   i  =  1,2,. ..,N!,   s  =  0,1,...,N,   ^  =  0,1,2. 

2k 


We   define 

N-s 


(15)  z(;;'^)  =  I  I  (A^Jt(-^))(^^^) 


i=0   j  =  6'^ 


®     _         V       7(^.3) 

,N  -       L     "^^-^ 

u-s>0 
or  even 


B=    (B^M'^',       B^^  =  0^(B^^) 
Q  =    (Q^V^)^/^       P  =    (pip^)l/^ 


4.2   Estimates  for  the   L„   Norms  of  the  Second  Derivatives 

Our  first  task  will  be  to  derive  estimates  for  the  Lp 
norms  of  the  first  and  second  derivatives  of  the  stresses  in 
terms  of  the  state  of  stress — t,  assumed  bounded,  and  the 
quantities  P,  Q,  and  Z.  We  will  not  at  this  time  delve 
into  the  structure  of  these  quantities  which  contain  terms 
due  to  the  non-linearity  of  the  constituent  equations,  the 
geometry  of  the  shell,  and  the  variation  in  thickness. 

Since  we  have  assumed  no  conditions  on  the  stresses  and 
their  derivatives  on  the  edges  of  the  region  M  given  by 
yttyO  _  1^  ^   ^g  introduce  a  set  of  auxiliary  functions  $., 
i  =  0,1,2,3,-.-   defined  by 

^0'^1'^2  =  (I-'^'^UV)^  =  (1.X-^U°^°)^ 
for  UV  =  U^J"  <  X^ 

K  —  — 

for  ]fV^   =  U^u"  <  X^4^"^,  where  k  =  5,4,., 
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0=0  otherwise  for  i  =  0,1,2 ,J>-  ■  -  - 

It  will  be  understood  that  t>     will  denote  any  particular 

t.      when  no  confusion  can  result.   It  is  easily  seen  that 
1 

1*1  <i,  i*;°'i  =i\\j  =ii*,j  '  &(v^>"')  =  6^(x"^). 

l*i^N^  '  =  '^i^N^  '  "  °  if   s  >  0,   i  =  0,1,2,..., N  =■■  1,2, 

For  any  vector  w,   defined  in  the  region  M  given  by 
U°U°  <  X^,   |U^|  <  H  we  define  the  Lg  norm 

llwll  =  Vj;||w|^dU-^dU^dU^. 
M 

LEMMA. 

(14)  lln(^^(°^l  =  ff(!|n°^  ,J    +  y  "^t^^    •'  +  >'^"t"). 

PROOF : 

From  the  appendix  (A2)  we  have 

(15)  Pl^\^^   =  0'("^i^',aa"  "■   X"^"i")- 
From  the  identities : 


'see  Sensenig  [6]. 
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^      ,aa  "  -     ,g2_  ■  -     ,g2_  "^  -     ,22' 

^     ,aa  ~  -     ,32        ^     ,21  ^  -^     ,§2' 

.22  _    ,23  ip  .33 

i      ,aa   -  ^      ,&2    ■   -      ,^       -      ,11' 

we  have 

(16)         i|*t'>'8_^ii  =  &(  I  iptj\^\  .  !!*t-_^||]). 

Setting  3=5  and  summing  (^.15),  substituting  in  (4,l6) 

gives  (^.15)- 

We  introduce  the  quantity  A  „  defined  by 

2        1 

,      .  _   .3k  .31c  a/2  _     Y      Yt(p)(s) 

Un    A^jj  -    it      ^ia,...a„T-     ,ia,...a      /  ~     L       L-        ,N 

± ^2Z±  i ilZi  u=l    s=0 


Denoting 


JJJ[.]dir^dU^dU^     by     J[']dU,      we  have 
M 


,,(1)(0)|,_     „^t(i'(|'ii,     i|«(2"o),|^     11^, 


t'^'!°'ii]  = 


and  clearly  we  have 
(18) 

From  equations  (4.6)  and  (4.7) 

^33         ^     ^3a       ^     _  3  ^03         ^  _..3a  a     _  a 

^     ,33       "^     ,3a+^,3^  ,3'     ^     ,37         ^     ,37       >7  "^  ,7' 

(19)   ti\33  =  -(l-v)(t5^^^^.t55^33)   +   (l+v)(Q3^-B5^), 

-  t^^ 
,aa  ,33* 


33         _  A+v\q     _„     ^   _  t^3         .  ^li 
*      ,aa  -  V^+V^^rr-Brr^        ^      ,r-        ^ 
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so  that  transforming  the  stresses  to  U  coordinates  and 
using  (4.13) 

(20)  ll$t^^  ,,"  =  OC'M  Jl  +  H{P^   ^+B^  ,)!'  +  Hz^^^n, 

(21)  ":/)t^^  ^^11  =  (J("0A^2"  +  "^^^^,33"  "^  "^(^+^^3)''  +  'i^Z^s"^' 

(22)  l'*t^^,aa"  =  G("^(Q+Brr+P^,3)''  +  "^A^^H  +  "^Z^g"^' 

Setting     h  =  a,      k  =  3      in  the   compatibility  condition 
(1.17)    and   solving   for     t  ^   ^,      gives 

(23)  Ht^  „''   =   CTC'^l''^   ^"+l'^tii      _II+IIM0_+B  „)"+II;^Z%t|). 


> '- 


,33  v/ V    ^^      ,^7      '^-      ,ag_       "^"^0^      ap 

y(  u )  ( s )  e 

t  ^         o      =   t    „ 

u-sX) 
or  even 

(24)  ll^ty    =   0("^(A,2   +  Z^2   ^   ^op  -^   ^rr)"- 
By   (a8)   of   the  appendix 

X"^ll0ht^°^^^^'l   =  D(x"^l!t'l   +  HX"^(l!0t^„n))    +  a(ll*t%l|). 
—        '_  ~"  — '_  '  — 

The  equilibrium  equation  (4.7)  gives 

Consequently  we  can  estimate  all  second  derivatives  of 
the  stresses 

(25)  X"^!l0ht^2"  =  (ji^'^M   +  Hx"^Ht^2'l  +  x"^"0hZ^2" 

+  x''^|l0h(B^^+Q)ll)  +  a("n^2''^' 
Using  the  equilibrium  equations  (4.6-4.7),  estimates  for  the 
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derivatives  of  h     (4.11)   and  (A3)  and  (A8)  we  have  for 

,01 

U-s  >  0  and  all  first  order  derivatives 

i\^t^\^)is)^i   =6'(x"^HI1  +  l'0(&fP)|l)  +  a(x!l0t^^^(°^!l), 
—    *i  —  —    *_ 

X"^l!0ht   II  =^(x"^l!t!!)  +  CT(Hllt%ll). 

,1.  ,  c. 


(26) 


It  remains  to  estimate  A  „.   We  make  use  of  the  fol- 
lowing  identity  for  T,Y,g  functions  of  Ijr,   twice  contin- 
uously differentiable,   g  a  function  of  ^f'     alone  and  sup- 
port of  g  contained  within  the  region  M. 


M  M  M 


Tg   T   .Y.-g.T   Y.+ 
J°,a  ,ai_  ,j^   ^,1^  ,aa  ,1_ 

M  ^ 


h  h 


Employing  (4.2?)  we  have  for   ll^A  ^ll^  from  (4.1?) 


M 


M 


-J 
M 


,a-  ,00-  ,i      ,a.-  ,a&-  >g. 


+0^^^ 


3k 


,a\ 


=  )   -  t 


3k 


,a 


+  t 


3k  /  *^ 


dU 


M 


X-^II^t^^f)  +  a('lM^2"^)> 


where  we  have  used  the  fact  that 
^(X"^),   and 


^,a^,a\l/2    /h  h  ^  \l/2 


,a  ,a 
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(29) 


Si  = 


,  aa    r^J  ,1> 


M 


.3k 


Substituting  for  Xr      ^.   in  (4.28)  from  the  equation  of 
compatibility  C^.S)  we  have 


(30) 


r^^t^^     t^^  ..du=  r^^t^^     j.du 

—  ,aa—  ,11—   J   —  ,aaK— 


M 


M 


1   ,^2,3k 


1+v 


^^^t^k   Uii     -  h  ut^i   .  IdU 
—  ,aaL—  ,3k    ,  o—   ,pkj  — 


M 


'   &(lJ^V\.aV2l  -  S. 


M 


+    ll-r^^     t^P  +  0^3k         Ttii        .  U     u\ii     IdUl), 

nL      ,S-      ,aa  -      ,aakj  ,3        \    ,P-/        ,P      -V 


M 


(31)      where      S,   =    V'^^'^X?'^        ^  ^^^  '-      -  ^^^^ 


o   -    1      ^    -  \t""   ,   -  Ih     u)t""      )\  ,dU|. 

2         'ji.      -      ,aa\        ,3        \    ,0-y        ,p/l,3  -' 


M 


Switching  the  order  of  differentiation  from  t_  , 

to   t   ,     and  integrating  by  parts  with  respect  to 

^        ,  KCtCt 

U        gives 

(32)  r*^^^  t-^  .,du  =  PrO'rn^^ .  ]  ^(t^^  ^  -'h  uVi   ) 

^-^    '      J      -      ,  gg—      ,^1   -       ^V-jL-      ,k2.        \         J 3a   I-    »£— '         » po/ 
M  M 


ii 


.ii 


^A2^3k 


+   21>        t^^   ,    -      h     u  it"'-'       )  dU|    f    1  ,0   t^  Jv^^l    +  S« 

,aV        ,3        i_    ,P-J        ,0/     -'         'J      -     ,ag  k  -'  2 


M 


M 
Using  the  equilibrium  equation  (^.7)  in  the  form 
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h 


(33) 


(3M 


3k 


,3 


h^^ 


—   ,ka    \  1^  —  /  ,a  ,a 


33 


,3 


•i'^^  ^  <^t'^)  ^  V, 


and  substituting  in  the  first  integral  appearing  on  the  right 

hand  side  of  the  estimate  (4.32)  we  get  for  that  integral 
h  „ 


(35) 


:-V2^t5n    +v      H(t^V^),    ^h  -  (h  u)tii 

M 

M 


ii 


h 


■,3 


3P 


M 

For  the  last  integral  in  (4.32)  we  have,  using  the  equilibrivun 
equations  in  the  form  of  (4.33-4.34),  that  it  is  bounded  by  an 
expression  of  the  form 
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(J6)   |jw^3t»_,„t"_^  .  (**,.)  ,6i'' 


M 


h 


+   0^Xh"^ 


h  . 


,aa' 


,  a\  h  /    ,  a  FT 


hx-^t^i 


,5 


■2''--   •  ^^^^  ""'-5  ^  ^^r^^^  v.ii  it^p 


-  .^-T^^    ,an'^^a-lI^ 


h 


h   h 


h-V^FT^  ,3a  ^^FT^  ,33  ^  "Tr:^  ,37]^^' 


+  S), ,   where 


Sh  =  Ilk*  «t 


'4   'J 

M 


.3P 


33 


g_—  ,aa     \\  h  —  /,a   —  ,3a/  h 


-  t 


33 


,a 


h 


^ -^-i'^, -^'*0,3^iii 


Substituting  the  estimates  (4.32),  (4.35),  (4-36)  into 
the  estimate  (^-30),  and  replacing  t    and  its  derivatives 
by  their  components  in  the  U  coordinate  system,  we  have  in 
terms  of  our  notation 


M  M 

+  Sg  +  s^  +  s^  +  x"^"n^i''^  +  x'^'itn^  +  iifi-V°^i'^) 


+  ai"n0t  ^'1^  +  Il0t^^ll^  +  x'^'lA^t  ^11^). 


-,2 


-,2 
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We  now  estimate  the  term  involving  J,   appearing  in 
(^•37)-   From  (4.8)  we  see  that  Jj^  contains  terms  due  to 
the  non-linearity  of  the  equations,  the  presence  of  body- 
forces,  the  variation  of  thickness  in  the  shell.   Transforming 
all  stress  tensors  and  their  derivatives  to  the   U  coordinate 
system,  and  using  the  identity  (4.27)  to  interchange  i_  and  a 
differentiations  we  get 

M 

—     '.   —  ,aa   ,2  -     1    5         ,£ 
M 

h  ^ 

where   S^  =  |  p^^f^t^^   t^H  ^^U. 
5    '^   L  h  -  ,aa-  J,3  - 


M 


We  have  for  the  integral  involving  the  body  force 

.h  _.         h 


+ 

M 


(39)   where  S^  =  I  j[t^^c.a^  ],  3^3111 


M 

Combining  the  results  of  estim.ates  (4,26)  through  (4.39) 


and  noting  that   A  =  [J  {^6^),      sup  A  =  U  (©),   we  have 

M 

+  I10B^J!^  +  il^J^B  J!^  +  l!$(Q+Xh"V°^!!^  +  X-^lltJ!^) 


..2 


+  a(x"^'|iht^2n^  +  Hty^^) 


33 


Hence,  from  the  estimates  (4.l8),  (4.24),  we  have 

(40)   lin^g"^  =  S^C'^TQ  +  Xh'^v|°^  +  V^^  +  z^2  +  >^""^v^!'i 
+  X"^'it!!  +  9'i0t^^'l^  +  "*B^  '1^  +  H^B     1!^ 


+  S^  +  Sg  +  S,  +  Si,  +  S^  +  Sg)  +  c(x"^n0ht^2n2). 


From  (4.40)  and  (4.25)  we  have  estimates  for  the   L^   norms 
of  all  possible  second  derivatives  of  the  stress  vector  t_ 
in  terms  of  the  L^  norms  of  the  boundary  tenns,  the  stresses 
and  the  quantities   Q,  V,  B  and  their  derivatives.   Substi- 
tuting into  estimates  (4.26)  we  have  similar  estimates  for 
all  the  first  derivatives  of  t. 


4.5   Estimates  for  Boundary  Terms  and  Body  Forces 
Occurring  in  Estimates  for  the  L^      Norms 
of  the  Derivatives  of  the  Stresses 

Boundary  terms  appearing  in  our  estimates  are  of  the 
general  form 


(41) 


J(fg)  ^dU  =  JJ[(fg)'^  -  (fg)"]dU^dU^, 


M 


u°u°<: 


where   f  =  f 


U^=+H 


and   f  =  f 


3 


We  use  the  inequality  for  positive  constants  K  and   D; 


,+  + 


lr/^+.r^-\,     +        -' 


,+  ^-\/  +.  -' 


(42)  |(fV-f"g")l  =  ^[(f%f  )(g  -g")  +  (f  -f")(g%g")] 
<  K(f^+f")^  +  K"^(g+-g-)^  +  D(f"*"-f)2  +  D-^g-^+g-)^ 
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In  our  estimate  the  quantity  S  ^     defined  by 

N 

■  '""    ■""  ■■'"  ^-Z  -1 


i=0 

+  ft(^)-^^^)+t^^)-(^)\'^X-^   +  rt(2)  +  (0)^j2)-(0)\-il 
+  \L        ,N-i  -    .N-iy"   ^   ^  V-    ,N-1  -    ,N-i/^   J' 


for  N  =  2,   will  play  a  central  role. 

From  (4.29)  we  have  using  (4.4l),  (4.42) 


1    'ov   —  ,aa—  »2/'2  ~ 
M 

irir<x 

+  l!J^'   (i''^gg-r^^gg)'""VVlH 
^\n^        (i^^^gg-r^^gg)''"'^IIVlH) 

^  e^-^!!.t^5^3|l^)) 


In  order  to  estimate  S^, ■ ' • ,S^      we  transform  all  stresses 
and  their  derivatives  appearing  in  these  estimates  to  their 
components  in  the   U  coordinate  system.   By  a  technique 
similar  to  that  used  for  estimating  S, ,   we  have  from  (4.^1) 
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Continuing  in  this  fashion,  it  is  easy  to  show  that 

+  >"^ll>ht  „"^  +  Ht^„ll^ 

where  we  have  estimated  all  the  first  order  derivatives 
appearing  on  the  right  hand  side  by  (4.26). 

In  order  to  estimate  the  terms  involving  the  body  forces, 
we  define 

N     -  Nl/2 


(^5)        B 


i=0    '   '-^l^'-^N-i   '^1 'N-i'^ 


(^b(u)(s).  b(,j)(s)^)  _ 


Prom  our  assumption  on  the  body  forces  given  in  terms  of 
(5-5)  and  the  estimate   a  =[J{^9^)      we  conclude  that 

U(ex"  )  if  u-s   is  even, 

gr(ex'^(9Qh/x))   if  u-s  >  0, 

fhB^jj_^;   hB^j^  =   Cr(e>"^^)  otherwise. 

Consequently, 

(47)    Ch^B^?^'!,  n^fiB   'I,  'I0B^?^'!,  I'^nB   I',  X"0h'^B^^^^?^'l, 

1|0b(i)|i)-,  n0B(^ni)|l)  =  9^(EX-'l!^lp. 

For  the  boundary  term  Sg  we  have  from  (4.39) 
S5  =  C^("*S^2"^  +  (ex'^)^"0"^). 
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^ .  h      Estimates  for  Terins  Arising  from  the  Transformation 
of  Lagrange  Coordinate  Systems 

For  the  expression  V  appearing  in  the  equilibrium 
equation  (^^.7)  we  have  by  (4.46) 

>||th-^v(?^"  =  &(>."i)h-^(B^°Up)^?^"  +  x"^!!^(&fP)i!  +  x-^i't") 

(^9)         =  OC'v^P^?^'^  +  x"^'!iP^'  +  X-^"ti'  +  (e/x^)ll0n), 

!1^V   ''  =  Er(x"^"iP!!  +  'I0P  J'  +  x'^l't'lU  a(l'0t^^"). 

fu  s) 

Investigating  the  terms   Z^liC'     (4.13)  which  absorb 

the  terms  induced  by  the  change  of  Lagrange  coordinate  systems, 
we  find  that  if   (u-s)   is  even,  the  general  term  in  Z  „  is 


of  the  form 


u^^)y-,^-'^^'y 


where   M^  <  N^,   M^  <  N^,   N^+N^  =  N,   M^+M^  =    Us .      If 

(ia-j)-M„   is  odd,  we  have   (i+j)-M,   is  odd.   For  if  not, 

adding,  we  have   ^i+i-M-,  -M^  =  ij-s,   v/hich  is  even,  a 

contradiction.   Similarly,  if   (n-j)-M„   is  even,  then 

(^j)-M,   is  even.   In  addition,  we  note  that  all  N-order 

derivatives  in  Z  ^  have  at  least  one  factor  of  A  multi- 
plying them.   Next,  we  observe  from.  (4.10)  that 
h^V    =  A^p^  =  0.   We  have  for  s  <  M,  N  >  0: 

(A^  %K    A^  %h    =   ff((nao)N-X-^)   otherwise. 
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Hence  if  N-s   is  odd,  we  have 

Consequently  if  Z^^'°'   is  such  that   (u-s)   is  even  or 
greater  than  zero,  then  every  term  in  the  sum  will  be 
maiorized  by  some  expression  of  the  form 

9  fix  -^t^^^  ^y   for  x-y  <  0  and  odd,    or 
u  ,1^ 

X  ^t}-^'    ^^'  for  x-y  >  0  or  even,  where  N^+Ng  =  N. 

In  addition,  all  N-order  derivatives  of  the  stresses   t_ 
appearing  in  Z  „  will  have  at  least  one  power  of  Qd 
multiplying  them. 

In  particular,  we  have  for  N  =  2 

(50)  ll^z^gii  =  ^{pne^^^W   +  x"^||n^il|  +  x'^H^ne^^^H  +  x"^l|tll) 


4.5  Analysis  of  the  Non-Linear  Terms  Occurring  in  the 
Estimates  for  the  Derivatives  of  the  Stresses 

Transforming  all  the  tensors   t  and  derivatives  to  the 
Lagrangian  coordinate  system  U   we  have: 

t  =  F(A)[t], 

f  =  t(J)  =F(A)rt|?)  +  ^[lh■^   At^l] 

t»  =  t   =  F(A)[t  ,  +  A  .t], 
»  *—  '  — 

where  F  stands  for  a  different  expression  in  each  case. 
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On  the  basis  of  our  assumptions  on  the  curvature  of  the 
undeformed  midsurface  we  have  from  (5.8) 

(51a)  Ti^i  =  ©(^oHV  ^),...,T>^jj  =  &(eoHVx^"^-^), 

,(0)"=  &(eo.^HVx' ),...,  7°^  =  &(eo^HVxN^^), 
jo)  ^  ej(eQQHVx^),...,.;S^  =  er(aonHVx^-^^). 

For  convenience  we  also  list:  " 

(5ib)(A|?).  .|,^)<A,O')=0(e,x-^), 

(r,(°);  ,(1^  .(^))=  0(eoHVxN^^),  l<s  <N, 
^  ,N    ,N    ,N  ^       ^ 

(riA^l)  =  &(9o^"*/^^){^,l\l^G,l)  '    '^^!l)  "  C5'(eQHVx^). 

These  relations  also  hold  if  all  derivatives  are  with  respect 
to  the  corresponding  U  coordinates.   This  is  easily- 
established  by  use  of  the  chain  rule. 

Transforming  P   and  Q    to  U   coordinate  system, 
and  using  the  estimates  (4.51)  we  have  from  (1.1?) 

pi  =  F(Ti,A,t)[tt  T  +  A  .t^  +  iTt  ,  +  nA  -,t] 

(52)    P^  =  F(Ti,A,t)[C^t^^  +  C2t],   where  C^  =  tT"  ( e^H^h/x^ ) 

and  Cg  =  &(0QH*h/x-^). 
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(53)  Q  =  F(r,A,t)rtt  2  +  ^  1^,1^  +  (^,1^^-^  "^  ^1-,1 

+  r  ^A  -|_t  +  (r  -l)  t_  +  (t  ^)   +  Ti  gi  +  T^t  2 

2 

+  rA  ,t  ,  +  A  „t   +  rA  ot1 . 

Differentiating  (4.52)  with  respect  to   U   we  have 
P^?^  =  F(T^,A,t)r  (t  +  T^)t  „  +  (A  -.t  +  t  ,  +  n  .  +  rA  Jt^?^ 

'i  —    —     '_      '  i.""     >—  '—  '—       '  — 

From  the  estimate  (4.51)  we  have 

Xh-^p|°^  =  F(r,A.t)[C^t^2  +  Cgtjj^  +  >^~\l^S?^ 
+  C^t  ^  +  C^t],   where   C^  =  ^{e^E^/y), 

Cp  =  &(1/X),   C^  =  6(9oHVx^)   and   Cj^  =  t7(l/x^) 
For  a  general  derivative  of   P  we  have 

(54)  P^^  =  F(r,A,t)[tt^2  +  A^^t^^t  +  (A^i)^t^  +  ^,ii,i 

2         2 
+    r   ^h   -^t   +    {r   -^)    t   +    {t    ^)      +  n  2^  +  i  2^ 

+  TiA  li  1  +  A  2-^  "^  ^^  2-^ 

=  F(r,A,t)[C^t  2  ■*"  "^2-  1  "^  ^-,1^^  ^  S-^'   ^'^^'"^ 

^1  =  D(Qo^HVx^),  Cg  =  6^(0oHVx^),  c^  =  6r(eQHVx^). 

We  define  for   N  >  2 

N-1  N-2 

(55)  0,„  =  I   X-1--P,,  .  I  ,'*'-\, 

i=0  1=0 
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Consequently,  we  have  from  (4.52),  (4,54)  and  the  fact 

that   P  ,  -  &(Q), 

1 
^-^^   ^,2  "   I  ^'""^^,1  ^  ^  "  F(r,A,t)[C^t^2  +  ^2-,  1  "^  ^-,1^^  ^  ^-^  ' 

where   C^  =  Q^(9QhHVx^),   C^  =  '^{Qq^*/1?),      and   C,  =  ©(e^H^/X^) 


We  investigate  the  non-linearity  occurring  in  the  boundary 
expression  B  „.   If  S"""   represents  the  surface  force  per 
unit  area  of  the  deformed  surface  in  the  Lagrange  coordinate 
system,  we  have  from  (2.4)  on  the  upper  and  lower  faces  of  the 
shell 

(57)         S^^  =  (^'^'T\a^'"')[^',3?r]^ 

where  f  are  the  direction  cosines  of  the  normal  to  the 
deformed  faces.  We  have  in  terms  of  £^ -ji,  the  skew  sym- 
metric e-system  tensor 

From  (1.2a) and  (4.58)  we  have 


=  KT 


By  the  tensor  transformation  rules  of  (4.1)  we  have 

-m3  ^  !^H,r3  .  U^  ^     ,ra\ 
77y\h      h  h  ,a   /^ 


3U^  -'       h  ^  '^   ^' 


or,  multiplying  by  -2^  |j,   and  summing 


au^h 


t3 


h-m3  au^  =  ,P3  .  ^   |_,pa_ 
H    aU^  ,a  H 
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On  the  surfaces  given  by  U  =  ±H  we  have,  using  the  notation 
of  (4.3) 


tP3±  _  (^^^±_   V,  ■H^±^  -  qP± 


^     +       ,a  ' 


=   S' 


-3\l/2 


g 


3/ 


=  (S^^  ±  h^pS^=^) 


'^\l/2 


g. 


From  (1.2)  and  (l.l4)  we  have 

i^g^J  =  1  +  F(n,t)[TTft]  =  1  +  F(in,A,t)[n4-tl. 

From  the  expression  for  S  „   (4.45)  and  our  assumptions 

»_ 

about  the  surface  forces  enumerated  in  (2.4)  we  have 


,2 


S  o  =  S*   +  Oiex'^), 
>  _ 

_3 
where  by  the  structure  of  g,  and  (2.4) 


S*,2  =  0"(C,2) 


4.6  Estimates  for  the 


\ 


Norms  of  the  First  and  Second 


Derivatives  of  the  Stress  Components  t^ 

We  begin  by  estimating  C  r>   in  terms  of  the   Iv>  norms 
of  the  derivatives  of  the  stresses.   We  have  from  (A4 )  and 
(A9)  of  the  appendix 

ll*i^t\I!   =  Cr(e[||</>t%|l   +   X-^||t|l]), 

ll^t   ^t     h|l   =  &(e[Hl|0t%ll   +    (H*)-^||t||]). 
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From  (^.56)  and  (4.6o)  we  have,  using  (3.6) 
From  (4.53a)  and  (4. 60)  we  have,  using  (3.6) 


+  x|l*h-H^^t^^|i  +  x"^|lt!|). 


^  Schwarz's   inequality  v/e  have 


Xll^h-H  ,t^,l!  =  S^(i'0t  ,t  J!)  +  a(X^H*-2!i0t^.t^  l|). 
Using  (4.60)  and  (3-6)  gives 

Xllrn-H^^t^J!  =  O(£(l|0t^2"  -"   H*"^''ill))  +  a(eX^H*)-2ll0t^2'!) 

=  &(ex'^ll*ht^2!'  +  ^'^llilD  +  a(ll0t^2")- 

Hence,  using  (4.26)  to  estimate  the  first  derivatives  of  t_, 

[xl!0h-^p^il!,  !!0c_2l1]  =  ?y(x"^ll0ht^2!l  +  ^'^":^" 

+  ^oll^zll)  -^  a(ll0t^2ll). 


To  estimate  t  ^  we  substitute  into  the  estimate  (4.4o), 
the  various  estimates  for  S^,...,S^   (4.44),   B  (4.45-4?), 
Sg  (4.46),  V  (4.49),  Z^2  (^-50),  S ^^     (4.59),  C ^^     and 
P®   (4. 61).   Next,  using  the  estimates  for  the  first  derivatives 
of  t   (4.26)  we  have 

il^t^ii  =  Cr(x"^i|tii  +  x-^!|0ht  11  +  ex-^li^l! 

^o"*i!2ll)  -^  a(ll0t^2")- 


t3 


To  estimate  t  „  we  substitute  into  the  estimate  (4.25) 
and  get 

Substituting  (4.63)  into  (4.62)  for  the  terra  x~''""0t  ^Jl  and 
choosing  0^  sufficiently  small,  we  have 

(X"^ll0ht^2''*"*i!2''^  =  &(>'^:'i"  +  eX'^'l^ll). 

From  (4.26)  we  have  for  the  norms  of  the  first  derivatives  of 
the  stresses 

(X'^l!0ht   'l,!!0t^  i|)  =  £r(x"^H!'  +  ex"^"^ll), 
substituting  into  (4.6o)  we  have 

ll*t^li%ll   =  &{ey-^{ep\\   +    lltli)), 
ll*h(t   j2|'   =    6^(eH*-^(llt|1+e!l0ll)),    HA  Jl   =    £r(ex"^ll*!!+x'^!ltli) . 
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Chapter  5  -  Estimates  for  the  L^   Norms  of  the  N-order 

Derivatives  of  the  Stresses 


Our  procedure  in  estimating  higher  order  derivatives  of 

the  stresses  hinges  on  the  fact  that  much  of  the  second  order 

ik 

analysis  can  be  adapted.   We  find  that  t   ^      for  fixed  6 

ik 

will  satisfy  basically  the  same  set  of  equations  as  the  t_ 

did,  provided  we  replace   P,Q,Z  ^V, . . .   by  P^?SQ  ?  ,Z  ,, 

V^ -,',....   The  equation  of  compatibility  can  be  used  to 
>  ±_ 

3 
estimate  any  derivative  of  even  order  with  respect  to  U 

1        2 

by  a  combination  of  derivatives  with  respect  to  U   and  U  . 

In  addition,  we  will  find  that  differentiation  with  respect 
to  U*^  of  the  stresses  or  their  derivatives  of  any  order 
introduces  a  factor  of   X    in  our  estimates. 

As  a  beginning,  we  recall  the  definitions  of  the 
quantities  t}^^  [l^      (^-9),    zj^J'^^   (^.12),  A^^     (4.1?), 
S^jj  (4.43),   B^jj  (4.45),  and  C ^^     (4.55)-   For  the 

quantity  t^^'^?fS  N  >  2  we  can  consider  three  cases 
(i)  la-s  <  0  and  odd,   (ii)  ^-s  >  0,   (iii)  ^-s  <  0  and 
even.   For  case  (i)   jj-s  <  0  and  odd,  we  have  from  the 
equilibrium  equations  in  the  U  coordinate  system  (l.l6), 
transforming  to  U  coordinate  system  and  introducing  Z  „ 

(1)      i-Hi>ht}^'>[lh\  =  Cy(x-^llW°)||-^^^l!  +  x-^l!*hc^^!l 

+  x"^ll0hZ^jj!J  +  X"^lI*hB^jj_i|!). 
By  (A8)  of  the  appendix  we  have,  provided   (s-^i)  <  N 


^5 


If  s-n  =  N  (>2)   we  employ  the  compatibility  equation  in  the 
same  way  as  was  done  in  (4.22-23)  to  reduce   (s-|i)   by  2. 
That  is,  for  general   r  >  2  we  have 

o)  x-^ii*ht(o)(^)ii  =  er(x-^!i0ht(o);^-2)ii  +  x-^ii0ht(2)(-2)„ 

+  X'-^ll*hC^jj|l  +  x"-^ll0hZ^jj!l  +  x"-^ll0hB^jj_^!l). 

Combining  the  estimates  (5.1-3)   we  have  for  u-s  <  0  and 
odd 

X-^||W^)(|)||  =  er(XH-l!|^t(0)(r-l),,  ^  |l^t^O)(r+l)„ 
+  X'^P^hZ^jjII  +  X"^'|t(°^[^:^^l!  +  X"^ll0C^jjll 
+  x"-'-"0hB  j^_^l|)  +  0  (ll0t^°^^Jj'-'-^!l),   for  some  r  odd. 
where  0<r<N,   N>2. 
For  case  (ii)  if   (n-s)  >  0, 

For  case  (iii)  if   (u-s)  <  0  and  even,  we  have  from  the 
equilibrium  equations 

r=l    '~  ' 
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By  the  compatibility  equation,  we  have  for  s-n  >  2 

Continuing  in  this  fashion  to  estimate  t ^  ' ^^'   as 
long  as  p  >  2,   we  eventually  arrive  at  the  estimate 

If  we  differentiate  the  equations  involved  in  the  estimates 
of  t^^^^2^   by  U  ■"■,..., U  ^'^   (see  4.14-23),  we  have 

+  ll*(z|S'°'  ^  z(|'0'  +  z(|'2)  +  z(|'2))|i 

Consequently, 

(5)  il*t^„||  =  e(ll*(C,jj  *   A^j,)||  +    ^  (  ^  IIOB^^Il 

^-sX)  r=l 
or  even 

(6)  X-^l^ht^jjll  =  S^(^o(1'<n"  -^  ^"^iIi!N^"^ 

+  x"^ll0h(C^j^  +  Z^j^+  B^jj_i)ll)  +  a(|l0t^^")- 

In  estimating  A  „  we  follow  the  procedure  used  in 
estimating  A  „.   The  most  important  ingredient  is  the 

47 


compatibility  equation  (4.8  )  and  the  equilibrium  equations 
(^03)«   Integrating  by  parts  numerous  times,  we  obtain: 

N 


i=l 

+  (sup  a)'10A^jj1I  +  !|0C^j^l|  +   "-^ll^ht^N'O 
^  N       ~        ~  ~    . 


where   sup  A  =   (0q)   by  (3.7). 


Substituting  the  estimate  for  S  „  and  for  the  body  forces 
from  i^.^j)   and  (4.46),  and  using  the  method  of  (4.49)  to 
construct  an  estimate  for  ^  „  .,   we  have  by  (5 '4)  and  (5-6), 

i=l 


Altogether,  we  have  from  (5*5  *  5-7  )  and  after  expressing 


Z  „  in  terms  of  the  stresses. 


N 


(8)  "nV  =  ^("*V  ^  (I  ^"'"^!ihi" 


+  «o'l*i!N"  '  +  ^^ 


•Vll)-o(Xx^-^!l*h-ipJg!,ll) 


i=l 


(9)     X'^jl^ht^jj"  ='-'(^oI^"^^,n"   +  I^":^!N-i"^ 


+  ex'^ll^ll)  +  a{Hty\). 


i=0 


48 


In  the  preceding  estimates  we  have  always  used  the 
function  0   without  regard  to  which  of  the  set   i^-}   we 
have  chosen.   Now  we  use  the  notation 


110,  f!i,  =  ;::^(0,f)^du^du^du^, 

Mi 

{(U^,U^U^)|U°'U°  <  X^IU^I  <  HJ,   i  =  0,1.2. 


Ml  =  ^ 


2 

{(U^.u2,U^)jU°u''  <  4-:^,  lu'l  <  H],   i 

—   —   —    —  —     ,.x-c     ■  — 


We   observe   that      1   =    \J  {1> .  )      in   region     M.      j    <  i, 


and 


hence   lyn^  =  '^{H.wW    ). 

We  also  employ  the  convention  that  if  no  index  appears 
on  the  norm  sign  the  region  of  integration  is  given  by  the 
function  1) .  ,      e.g.,   i|0.|!  =  110 . 'I .  ,   and  for  example  from  (5.9) 

i=0 
.  +  eX-^^1101^1!)  +  a('l0t^^il). 

It  now  remains  to  estimate  the  terms,  such  as   C  ^, 
containing  the  non-linear  behavior  of  the  equations  of  equi- 
librium and  compatibility,  and  the  boundary  conditions. 

Generally,  we  have  that   P  „  and  Q  ^  -i   have  the  same 

form  as  an  expression  in  the  quantities   (A,Ti,t_)   and  their 

derivatives.   The  same  is  true  for   X~  p    and   P      and 

,s_        ,  r+s 

similarly  for   X  Q    and   Q  „   .   From  (4.55)  we  calculate 
C  ,  using  (4.52)  and  (4.53).   We  find  that 

C^,  =  F(A,Ti,t)[C^t^^  +  t^it^3  +  (l,i)"  +  ^2^^,1^^ 
+  C,t_  2  +  C^t_  ^   +  C(-t_],   where 
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Noting  that   H  <  H(H*  )""'-h(U-'-,u'^ )  ,   we  have 

(10)    Hi;0.C  ,'!  =  Q'Ce^^l'^^ht  ,''  +  H(H*)-^("0,ht   t  ^" 

+   "0,h(t^^)-"  +  X'^"*3h(t^^)2|!) 
+  9QX~^;i0,ht  2'!  +  Qq\'^"^^^1   1"  +  9x"'l!^ht") 

It  is  necessnry  in  estimatinr^   C  ,   to  employ  a  notation 
for  the  pointwise  no'^s  of  the  derivatives  of  the  stresses. 
We  then  make  use  of  a  Sobolev  type  inequality  in  estimating 

the  non-linear  terms.   For  this  reason  v/e  introduce   T,   and 

(0) 
T^    as  follows  : 

(11        T,  =  max|0,t  ,h|,   T^    =  max |0,t ^ , ^ | , 


and  for  the  N-order  derivative 


(lib)  T^  =  max    l*N+2t   Nh|,    T^^^^    =  max    10^+2^1°^  I- 

^N+2  -  %+2 

Employing  this  notation,  we  write  using  the  estimate  (^.69) 

and  the  inequality  (3.6), 

(12)  A!l*jht_^t^2|l  <  T^^llht^^llj  <  T^^!'*,ht^3i; 

=  ^(W-/^\^),      where   |!ll|^  -  gTivTiX). 
Using  Sobolev' s  lemma  (AlO)  of  the  appendix,  gives 

(13)  T^HX"^  =  Q^dl^^t^^hl'  +  X'^i'tllo). 
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Similarly,  we  have 

Estimating  the  remaining  terms  in  (5. 10)  in  a  similar  fashion 
gives 

(15)  Hii0,c^^n  =  ©(e'^^ht^^il  +  e^x-^M^+   a(i!0,t^^|l). 
To  estimate  the   L^   norms 

!'0^h"^-p(°)li  and   |!0^h"^p(°^ll 

we  calculate  the  form  of  the  expression  P  o  •   From  (4.52) 
we  have 

(16)  p(|'  =  F(.,A,t)Lt(|'t_3  +  t^^t(|)  +  [t^^'f  t_^ 

+  Cf-t  ,  +  C^t}^'    +   C^t  ^  +  Cgt],   where 

c^  =  "^{e^E^/x^),     c^  =  &(i/x),  c^  =  &(eQH*h/x^), 
Cj,  =  &(0oHVx^),  c^  =  &(eoH*h/x),  Cg  =  &(0qHVx^), 
QH*h/x^),   Cg  =  &{eo^ 


C^  =  &(0^H*h/X^),   Co  =  &{enHVX^). 


We  note  from  (4.54a)  that  this  is  also  the  form  of   >~  P^^^ 
Consequently,  using  (5.11-5.14)  to  estimate  the  terms  in 
(5.16)  gives 
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+  x'^'l'^jtnj. 


Using  the  previous  results  for  N  <  2   (4.64-63)  we  have 

(17)  ll03h-lp(|)!!  -K  X-^ll^jh-lp;^),,  =0(Ti°^-^(/HX) 

+  T^X"'^\/HX)  +  ex"^(v1iX)  +  X"^l!0^ht^^"J. 
Using  the  formula  (AlO)  from  the  appendix  we  have 

(18)  t[0\/hx  =  ti(x^ll^3t||^!l  +  Hll^jht^y,!  +  X-^!|tllQ). 
Substituting  (5-13)  and  (5.l8)  into  (5.17)  gives 

(19)   "V--'p||^!l  +  x-^n*3h-ip(o)„  =  er(x-^!!03t(|)!| 

+  X'^l'^jht^^ll  +  x'^lltllQ  +  ex'^(vTTx)). 

From  (5.8)  we  have,  upon  substituting  the  estimates  (5"15)> 
(5.19),  (^.64),  (4.65), 

(20a)       H^\^^^\    =  ^(x'^ll^^ht^^n  +  ex"^(yHX) 
and  from  (5-9) 
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(20b)    X"^!*,ht^^|l  =9  (%[x"^"*3ht^3"  +  "^^i!^"] 

+  eX'^(v^X)  +  X'^llf'  )  +  a(!l0^t''^"). 
If  b        is  sufficiently  small,  we  have  from  (5-20) 
(21a)  ;>-.t^3!l  =  6(£X'^(yHx)  +  x'^iitllQ^  =  tJ'(ex"A/Hx), 
(21b)  X"^:>-.ht^^;|  =  0^(ex"^(yHX)  +  x'^iltH^)  =  &(ex"'yHx). 

Substit'.;:.1npj  the  estimates  (5-21)  into  (5-13)  snd  (5.l8)  we 
have 

U^ht^^i  <  T^  =  £?(£), 
|03t|^^'<T|°)  =  ^(sx-^). 

Continuinp:  this  same  process  by  estimating   C  .,   P  ^s  i 

>  i_    » -J-~ -*- 

anri  us  in";  a  function  0.   with  support  contained  inside  the 


repion  M.  ,,   v/e  will  get  pointwise  estimates  on  T.  _„   and 
T:^i, ,   (i  =  ^^5,...)'   Hence  by  choosing   Gq  sufficiently 
small  we  will  be  able  to  get  pointwise  estimates  for  as  high 


0 


a  de>\' "-^tive  of  the  stresses  as  we  require.   At  the  point   P 

1    2 

on  the  mid  surface  v;e  have   U  =  U   =  0  and   $.  =  1. 

Consequ'^mtly  8t   P^^ 

—        ,  U.-|  .  .  .  U.„      ^ 

for  as  high  a  value  of   N  as  mq   q.b.tt^   out  the  estimates. 

By  the  transformation  formulas  (4.3,4.5)  and  the  estimate 
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/h  h   =  g (0-(h/X)   we  have  from  (5-22) 


For  the  non-linear  terms,   P,Q  and  their  derivatives, 
we  have 

P  T  =  0(eX"^),   P  o  =  S'CeX"^),..., 
Q  =  0(£X-2),   Q  .  =&  (£X-5),.... 


5^ 


Chapter  6  -  Summary  of  the  Estimates  for  the  Derivatives 

of  the  Stresses  in  a  Thin  Shell 


From  the  differential  equations  (I.16-I7)  we  have  by  an 
analysis  similar  to  that  carried  out  in  (^.17-^.24) 


ex   ' 


Ij-s  >  0  or  even. 


^(eX  ^^^-"-h   )   M-s  <  0  and  odd. 


Bv  employing  the  boundary  conditions  we  can  improve  the 
estimates  on  certain  of  the  stresses  and  their  derivatives 
Since, 

t^^(U^)  =  t^"-  +  Jt^'^^dU^,   k  =  1,2,3 


•h 


we  have 


t^5  =  tr(Ba),...,tP' 


,a. 


■  a 


=  menx-^). 


N 


t'J  -  ^(t5'+ .  t"-)  =  e'(ea^). 


2,-N> 


, a^. . .a^ 


=  m^^">^")* 


for  as  high  an  N  as  we  have  carried  out  our  pointwise 
estimates. 

V/e  summarize  our  result  by 


(1)   t(")(s)  .  5^1(t(^)++t(2'-)(^)=  \ 


0(ex"^^)  ^.-s<_0  and 

even, 

0^(eX-%^-^)   u-s>0, 

&(sX-N-^VM  U-s<0, 
and  odd 


where  p  =  0,1,2;   s  =  0,1,.. .,N. 
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These  estimates  hold  at  a  point   P  of  the  shell  vmich  lies 
on  the  normal  to  the  unstrained  midsurface   2   at  a  point 
P  ,   in  the  parametric  coordinate  system  U   which  has  been 
identified  v/ith  the  unstrained  Cartesian  coordinate  system 
X  ,   which  has  its  orir'in  at   P^  and   X   axis  normal  to  2, 
at   Pq. 

V7e  can  -generalize  the  estimates  {6.1]    to  the  covariant  end 
mixed  components  of  t   by  noting  that 

t^^  =   t^^  +  F(-„t)[r.t]  =  tj  +  F(n,t)[Tit], 

and  r.^^,  =  ^{e^ny-^'^). 

Consequently,  it  makes  no  difference  in  the  estimate  (6.1) 
whether  the  stresses  are  (~i"^'en  in  contrava riant ,  covariant, 
or  mixed  components. 

By  using  (3.9)  and  an  inductive  argument  on  t^ie  order 
of  differentirtion,  we  can  sliov:  that 

Estimates  valid  in  a  particular  coordinate  system  are 
valid  in  all  coordinate  svstems  once  they  are  written  in  an 
invariant  notation.   '/e  shall  always  assume  that  whatever 
coordinate  system  is  used,   U*^  will  be  constant  along  normals 
(here  r;±ven   by  -h  ^  U  <^  h)   to  the  undeformed  surface. 
Under  these  conditions,  any  coordinate  system  U   can  be  used 
to  specify  estimates.   For  instance,  if  w.,   denotes  a 

covariant  tensor,  then  w  behaves  as  a  tensor  under  any  trans- 

1   2 
formation  from  one   (U  ,U  )   system,  to  another,   w,   behaves 
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as  a  vector,  and  w,-,  behaves  as  a  scalar.   For  the  stresses 

we  have  {t°^t^)^^^   =    f]{e) ,      (t^^t^")^/^  =  e^  (efi) , .  .  .  .   In 

1   "^ 
addition,  we  see  that  any  coordinate  system  U  ,U"   can  be 

picked  to  specify  the  conditions  on  the  body  and  surface  forces, 

Finally,  we  remark  that  the  Kirchhoff  hypothesis  that 

"normals  to  the  middle  surface  in  the  unstrained  state  go 

over  into  curves  normal  to  the  deformed  middle  surface"  is 

approximately  true  in  the  sense  that  the  angles  made  by  the 

a  3 

lines  dU  =  0  with  the  surface  dU  =  0  are  at  most  of 

order  eQ,      rather  than  of  order   e.   This  is  seen  by  noting 

that  the  cosine  of  this  angle  is  given  by 

^  3?/2   ,     3  3 


=  1  -  2(l+v)t^t^  +  t^(e^) 

=  1  +  er(eV). 
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Chapter  7  -  Description  of  the  Stress  System  and 
Geometry  of  the  Deformed  Shell 


We  choose  as  the  functions  for  the  approximate  description 
of  the  stress  system  in  the  shell  the  symmetric  quantities 

■z. 

where  we  write  u   for  U"^.   By  means  of  the  estimate  (6.1) 
we  have  that 

This  expresses  the  fact  that  along  a  normal  fiber  of  the 
shell  the  longitudinal  stress  varies  linearly  and  is  deter- 
mined by  T  ^   and   -  ^   within  an  order  of  error  of  magnitude 
•^   ogi        op 

eQ'^      at  most,   t  „   can  be  associated  with  the  resultant  stress 

ap 

on  a  fiber  coinciding  with  a  normal  to  2:   in  the  unstrained 

state,  and   o  q   with  the  moment  in  a  similar  fiber: 
op 

%^6u   =  2hT^^  +  9^(ehfi2), 

-h 

Jt^^udu  =  pa^p  +  9^(eh2Q2^. 
-h 

To  describe  the  geometry  of  the  deformed  shell,  we  use 

the  coefficients   e  ^   and   i  „   of  the  first  and  second 

dp        Op 

fundamental  forms  of  the  deformed  midsurface  referred  to  the 

1   2 
parametric  system   (U  ,U  ). 

(2)  ^a&   =    (6a6)u=0 
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-l/2_'  \     rr-JT^/^/      3 


(?'  V  =  ([^Sr  'r^)^.o  °  LtsjT  (i-ae  +  c^) 


u=0 


We  refer  all  tensor  operations  performed  on  the  surface 
2   to  the  metric  defined  by 

dS^  =  E  „du"dU^. 

By  (2.1)  we  have  on  the  surface  S  • 

g""^  =  G„^  =  0,   G^^  =  G^,  =  1. 

On  this  surface,  passage  from  the  covariant  to  contravariant 
components  is  the  same  as  that  defined  for  a  two-dimensional 

tensor.   For  example,  the  tensor  w,    gives  rise  to  the  two- 

a    3    a    a  _ 

dimensional  tensors  Wq„  ,  Wq„.  w-,„.  Wq,,...,   on  the  surface  ^^. 

p7   p7  ^y       P^  ^ 

Covariant  differentiation  with  respect  to  the  metric 
dS^  =  E^dU^dU^  will  be  defined  by  the  symbol   "1"  and  we 
will  use  the  Christoffel  symbols 

Covariant  differentiation  in  the  ds   =  e  ^du'^dlT  metric 
will  be  defined  by  the  symbol   "1|"  and  will  employ  the 
Christoffel  symbols 

a  — aa  ° 

where     Co        is   a  tensor,    and      e   ^e  „   =   5o. 
PY  ufj  p 

We  choose  to  derive  all  our  estimates  in  the  special 

1   ? 
U,U  -system  employed  before,  but  all  results  of  the  form 
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w^    5  ~  ^  ^^^   ^°^  ^■'•■^  components  of  the  tensor  can  be 
interpreted  to  stand  for  the  equally  valid  and  invariant 
statement  w  „    ^\i°^'^"'^   =  ^(Y^).   Thus  for  the  deriv- 
atives of  the  stresses,  we  have  that 

"h--   ^.Ti-.-Y^'S^'^^"")   ="^   ho^,7i...7,  -©~<^^''''' 

and  we  assume  these  statements  to  stand  for  the  analogous 

invariant  statements. 

The  quantities  i^      and   e^   (7.2-3)  will  describe  the 

deformed  midsurface  modulo  a  rigid  motion.   To  obtain  the 

displacement  of  the  original  undeformed  midsurface  2    we 

integrate  a  system  of  ordinary  differential  equations  due  to 

Gauss  for  x     and   e  . 
,a 

We  show  that 

By  rotating  and  translating  the  deformed  shell  we  can 
arrange  that  at   P- 

From  (1.1^)  we  have  for  any  point  in  a  neighborhood  of   P 


_ 

See  John  [2]  upon  which  this  analysis  is  based, 
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(det  x^^j^)^  =  det  Sii  =  1  +  &(£+fi^). 

Hence,   (det  g^^j^)"   will  be  bounded  for  small   0^. 

Conseauently,   (x^  ^  )"''",   defined  by   (x^  ■)~'^x^   v  =  ^i 
will  be  bounded.   Since 

we  can  now  conclude  that  x^  .  .  =  \j    ((e+fi  )/h). 
For  any  point  in  the  shell  we  have 

From  the  identity 

conclude  that   (x  t;  x  ^;  x-^  , )  =  1  +  (7(Q  ) . 


we 


From  the  identity 

((^^,i^'^  -  4^  +  (^\k  -  \)  =  -((^^i^"^  -  ^jH^^k  -  ^k)' 

and  the  fact  that   (x<^  .)"■'■( x-^'  .  )  =  5.   and   det  x"^'   =  1  +  &{Q^) 
we  have 

(x^^,)-i  =  5J  +  0(n),  {{x\^)-^,   {x\^)-\   {x\^)-^)   =  1  +  e^(fi^) 

From  (7.4)  and  (7.5)  we  have 
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However,   c^,  =  U(^  /x)   and  on  the  surface  u  =  ±h  we  have 
from  (4.57) 

f^  =  1  +  0(q2),   §^  =  &(n),   but 

5^  =  G^-^E.  =  5^?.  +  &(fi^),   and  consequently 

?^  =  1  +  ft^c.a^),  f  =  er(.'^). 

Setting  1=5  In  (7.7)  and  substituting  for  z       we  arrive 
at  (7.6). 

We  Introduce  the  quantity  o)^,   given  by 

2 

If  fi   is  of  the  same  order  of  magnitude  as   e,   we  shall 

have  occasion  to  introduce  the  quantity 
(9)  o)  =  x^  -  f. 


where  a)|^  =00^  +  ^{^^  {e+^^)/h) 
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Chapter  8  -  Interior  Differential  Equations  for 
Shells  of  Variable  Thickness  with 
Surface  and  Body  Forces 


We  define  pseudo  stresses  T.   by 
(1)         1?=  (6/G)^/2t%  -  W6^=  (W3_^-  W)6' 


+  (2W   +  2W   )e^  +  (4W   -  3W  )e\l 
^1  ^2   ^      ^2     s^   s  1 


-  W5^  +  F(t)t'^, 

where  W  =  4(l+v)t^t^  -  vt^t^  +  F(t)t^. 
?  s  r     r  s 

Using  the  stress  strain  relations  and  the  equations  of 
equilibrium  we  have 

5  ni        tnr  s        ms        l.ms     "1 
^rm^si  "  ^  ^rm^si      ^si;m  "  2^  ^sra;ij' 


^rs_ni 


s 
It  follows  from  the  definition  of  c    that 

rm 

We  put   down  for  reference, 

(M       (g/G)^/2    =   1   +    (l-2v)tj   +   F(t)t2 

gga  =   &sa(l+F(^''^)   LTT+t])    +  2(l+v)tg^  +  F(  m,  t )  [t^+r,t] 

6^ 


By  (8.4)  we  have  from  (8.3)  on  the  surface  u  =  0, 

Where   P.  =  2(h/x),   and   (\)u=o  "  \* 
From  (8.1)  we  have  on  the  surface  u  =  0, 
(6a)  T^  =  ^i  +  '^(e^),   and  by  (8.2) 

(6c)         "^a  =  "^a  ""  Q'^"^^^)  =  ?r(ea). 
We  have  for  some  intermediate  value  v 

Evaluating  this  expression  for  u  =  +h  and  u  =  -h  we  have 

^(^'  -  ^''  '  ('^.j'u.o  ^  K(^,353'u=V  - 

(<3'u»0  =  ^(Tr-^"'^^("''/^'- 
We  obtain  from  (8.5) 

Writing  the  third  equilibrium  equation  in  terms  of  the 
pseudo  stresses  we  have  for  u  =  0: 

(8)    %  -   ^^l,,K-0   =  ^3|a  -  ^^a  -  ^K  -   ^,3  ^  S 

=  ^^(n^en/x). 

Here  by  (8.6c)  and  (6.1) 
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Consequently,  we  have  using  a  Taylor  expansion 
(8a)  T^^+T^'"  fT^+4.T^-) 


Now  performing  covariant  differentiation  with  respect  to  the 
metric  in  the  undeformed  surface  gives 

+  gTlQ^ea/x),   where   (B^  .    ^  )„-o  =  ^^  i    i 

Expanding  the  quantity   (T^+W)   along  a  fiber  perpen- 
dicular to  2    we  have 

(9)    ^(Tf  .   w^)  -  (T^-  .  w)]  =  (i^^3  +  w^3)^^, 
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Calculating  T^  +  W  in  terms  of  the  stresses  t .  we  have 
from  (8.1) 

T^  +  W  =  t^  +  2(l+v)tgt^  +  (l-4v)tjt^  +  F(t)t^, 
using  an  identity  similar  to  (4.42)  we  have  by  (6.1) 

=  (tf-t^-)  ^  (i-^v)(t^^t-)[|(tf^t^-)]  +  f^ici^ea^). 

Since 

we  have 

(10)     ^(l|^  +  W+)  -  (T^-  +  W-)]  .  i(tf  -  tf) 

+  (i-''v)(t°_5)„J|(t^''  +  t^")]  +  e'(n^enVh). 

We  now  seek  an  expansion  for  (W  -z).._p.     in  terms  of  the 

1/2 
Stresses  and  the  boundary  conditions.   Dividing  by   (g/G) 

and  differentiating  (1.6)  with  respect  to  U^,   we  have 

Prom  the  expression  for   i  „   (7.3),  we  find  using  (8.4), 
(l.l4),  and  the  estimates  for  the  stresses  and  their 
derivatives  from  (6.1), 
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On  the  other  hand,  from  (l.l4,  6.1,  7.1,  7.2)  we  find 

Therefore,  from  (1.4)  we  have 

(12a)       0)^  =  -(l+v)o^  +  va^J^E^  +  ©'(J^^^/h), 

(I2h)  03^  =  .(i-v)a„„  +  ©"(n^eA). 


oa    ^    '^aa 


Prom  (8.6a) 

Upon  substituting  into  (8.11),  we  have 
Substituting  for  (t°  ,)„_n  =  ot     in  (8.10)  we  have  from  (8.12b) 

CX  J  ^   Cl  — V-*       t* 

W         ^(t|+  +  w+)  -  (^-  +  w-)]  =  ^(^^  -  t|-) 

Next  we  note  that  for  u  =  0: 
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=  -^%a;P33-'^,J3-^("'^/'^^'' 

Substituting  (8.15),  (8.14),  (8.13)  into  (8.9)  and  substi- 

3 
tuting  for  T^  ^  in  (8.8)  we  conclude  that 

The  Codazzi  equations  for  the  deformed  midsurface  and 
undeformed  midsurface  of  the  shell  are  respectively: 

We  have  in  terms  of  00  _   and   c  ^  : 
This  gives  us  two  non-trivial  equations : 

(")  <o„^l^  -  co^^ie  =  e?("'^Ax). 

We  write  the  Gauss  eauations  for  the  curvature  as 
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'^11^22  "  ■^12'^12  ~  2^^12,21  "*■  ^21,12  "  ®11,22  "  ^22,11^ 

"  2^^12121  "^  ^21|12  "  ®11|22  "  ®22|ll^ 

+  i^a^^lA2  -  ^12^21^  +  ®X^^''l2^21  "  ^11^22^' 

and  express  the  quantities   e^  in  terms  of  t^  and  Qj-q. 
For  this  purpose  we  employ  a  set  of  approximate  "constituent" 
equations.   From  (I.l4)  using  (6.1,  7.1,  7.2)  we  have 

Combining  the  constituent  equations  and  the  Gauss 
equations,  we  conclude  using  (8.7)  that 


^11^^22   +   ^-"^ 


22^11   "    ■4.2"^21 


-   L^.O),-   =   -(1+v)t 


■^21   12 


XX|ti|i 


XXIuiJ        ^2^    5  ^    ' 


UM 


-    (1-^v)[b^ 


3+     3- 
_^  t^  -t 

2h 


]l-l( 


'"ll"^22   -  ^12"^21 


)    +    ^'(fi^e/x^). 


Writing  the  equation  invariantly,  we  have 


/■,o\        X,u         a  ^       ^3  a       1,   a  p         a  p.        v/^3+       4.5\hM 


+   (1+v) 


B„  + 
a 


a       a 


a 


=  ^(n^eA^). 


"5+ 
It  remains  to  express  the  terms  tri   in  terms  of  the 

applied  surface  forces,  and  the  U-^  derivative  of  the  stresses 

,  a  ,       a  01 

(tc  -z)  r\  =  Oa     in  terms  of  the  curvature  changes  o)-, 
p,p'u=o    p  -°     p 

Prom  equations  (8.12)  and  (8.17)  we  have 
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(19) 
(20) 


^a& 


^  A 


r^rv  -  tiXyV  "" 


0  (^^^e/h). 


ay 


1  ,aT 


L+v 


CO 


y  -  ^"uh^""  *    ^(^'-/^^) 


1-v  ^ 


From  equation  (^.57)  we  have  for  the  surface  forces 
(21)   S*  =  t^*-  h^^tp-  +  ^(a^eQ/x), 

S^  -  S"  =  (t^""  -  t^")  -  h^^(t°'^  +  t^")  +  &  (.^^e.a^). 


Consequently, 


(22a) 


(t^^t^-) 
^  g   g 


3+  .3- 


(S+-S") 
2      h      2h 


(t  +t   )  h  Q    v-^^"^  )        r»      o 


(s!  -  S-) 


=  ^'a¥-^^V^^0^^^^^-/-)- 


Similarly,  we  have 

(22b)   |(t^''  +  t^")  =  |-(S^  +  S-)  +  6^(0^6), 

(22c)   ^(tp""  -  t°-)  =  ap  +  &{^^B/h), 

(22d)   ^(t^""  -  t^-)  =  ^(S^  -  s;)  +  o'^h^^h^p  +  er(a2e.Q2/h), 

^22e)  |(t^^  +  t^")  =  hh^po^  +  1(3^  +  s;)  +  ef(a2eQ). 

Inserting  equations  (8.19-22)  in  equation  (8.I6)  we  have, 
upon  rearranging  the  terms,  reintroducing  E,   Young's  modulus, 
and  expressing  the  body  forces  as  integrals 
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1-v 


.Cl/.P      ,         Px  1/0  + 


.°^  -    ^izlvv  «> 


+   T"(r   +  03^)    -   i(s;   +   S:)(L     -    (iZ^V)^-) 
p^   a         a'        2^   5  3       ot        ^  1-v  '   o.' 


+  7(s^  +  s;)r +  -t,(s 


;)i°-^(s; 


-S'  ^  arJ.^V" -^  ^^V^a"^""" 


t^( 


Ef^  eh/x    ). 


From  equations  (8.7)  and  (8.l8)  we  have  respectively 
(25b)  (hT^)lp  +  hB^  +  i(S^  -  S")  =  gf(En2eh/x),   a  =  1,2; 


We  complete  the  set  of  six  approximate  equations  for  the 

P        a 
unknowns  cd   and   t   with  the  equation 


(23d) 


P 


"al7-"7la=  S^("'^Ax), 


which  provides   two  non-trivial   equations. 

As   a   consequence  of    (8.8a)    and    (8.20)   we  have 


a, 


2 
h^  Y 


2(l-v^) 


^p   ^^p      ^  hi- 
7  2  +     2^P,3   + 


6r(Ea2eh/x) 


Since  t^  =  &(EeQ)   and  t^^^^^  =  B^(Ee/>^)   we  find  for  t 
the  approximation 


71 


2      2' 


y 


(sj  -  si)      (s+  +  S-) 


t^    - 


2      2\^ 


3 
For     t-,     we  have. 


tl  =  |(s^  +  S-)  +  g^CEa^e) 
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Chapter  9  -  Interior  Differential  Equations  for  Shallow  Shells 

of  Variable  Thickness  with  Restricted 
Body  and  Surf?^^  Forces 


We  shall  assume  that  there  exists  a  potential  function 

1   2 
S  =  S(U  ,U  ),   such  that  the  resultant  of  the  surface  trac- 
tions on  the  midsurface  is  given  by 

where, 

S  rt    «  =&  (e^/X^'"^).   P  =  0,...,N+1. 
For  the  other  surface  forces  we  have  from  (3-5) 


(S^+  -  S^-)  ^    ^   =&(£nVxP),   p  =  0,...,N 
'  1     p 

(S^""  +  S^-)  ^    ^  =&(e/xP),   p  =  0,...,N, 
'    1"  '    p 


P 

and  for  the  body  forces  perpendicular  to  the  midsurface  2 
we  have  from  (3.5),   B,  =  [f  {eQ/\)      and 

Q'(e/X^"*''^)   if  an  odd  number  of  the 


0 


^,il...i 


=  < 

P 


i  , . . . ,i   equals   3, 


^(e/x^h)   otherwise. 


We  assume  that  the  body  forces  in  the  plane  of  the  shell  are 
negligible 

B^  ^     1   =  ?r(e«VxP'^-'-),   p  =  0,...,N. 
a,i-L...ip 

From  (8.23)  we  have  under  the  assumptions  of  this  chapter. 
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for   Q  =  2(h/x) 

(la)   (hT^)lp  +  S^^  =  tKsa^eh/x), 


+  L  'Uq  -  L  ai_ 
a  3    a  P 


S"-  a  p    a  e 


(Id) 


-E   1 


h'.°)|^-(i.v)((h3,i^.^(h')i°.p); 


^|(sX)r 


If  we  assume  that  H/R  is  close  to  e  in  magnitude,  we 
can  derive  equations  for  shallow  shells  which  resemble  those 

Q 

of  V.  Karman  and  F5ppl  derived  for  plates  of  constant 
thickness. 

If  we  assume  an  Airy  stress  function  to  exist,  then  it 

Q 

will  satisfy  equation  (8.4a)  to  the  required  order  of  error. -^ 
For  <I>  the  Airy  stress  function  we  set 

,  ci    ,  a  Y    a  7,  ,  6     a 


3 


3  5    5  3'  '7 


^3' 


Performing  the  covariant  surface  differentiation,  we  have 

(^-^1)13  ■^^l  =  ^'21  -  ^ti2  =  ^1X212 

(hx^)lp   -.3^2   =   *ll2    -   ^'21   =  *lXl2r       • 

n 

See  V.    Karman   [4  ]. 

^See   Green   [1]. 
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where  R  qc-„   is  the  Riemann  tensor  for  the  surface  2^. 
apo'Y  0 


We  show  that 


^J'TR-.^c-.  =  S^iEa^CeQ+Q^)). 


a, 
'0357 


At   Pq,   we  have  ^^(Pq)  =  i^   a^^O^  ^  '^'   ^^^  ^^  ^^* 
^(Pq)  =  ^   ct(Po^  "^  "^^   -^  ^^^  assumptions  made  on   f   (2.3), 
we  have  ^Q^y   =  Cr(^'^^/X").   Hence,  it  will  be  sufficient  to 
show  th^-t 

(2)  .  <t!^  =  *^^  =  0(E£P.X^  +  Eil^X^). 

Expanding  <J    about  the  point   P„,   v-ze  have  for  some  inter- 
mediate  value   P' 

and  for  <t>  ^^     we  have 

Since  *la3(P)  =  S'iEeh),   we  have  <I'^ap(Po)  =  'S'CEeh).   The 
expansion 

together  with  the  fact  that  <I>  ^^  is  bounded,  gives 

ct>       =    CP(Eea  +  EX^). 
Substituting  this  result  into  (9.^)  gives,  using  (5.7t>)  to 


estimate   iJa}* 
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and   consequently  from   (9.3) 


^  ^  =    CTCEehx  +  Ex^a^). 
>  CI 

Substituting  into    (9-^)   gives 

*,ap   =    9(Eeh  +   EP/x), 

inserting  this  into  (9-3)  for  0     gives  (9-2). 

Recalling  from  (7.9)  that  cuj^p  "  ^op  "*■  ^"^^  /^^  " 
-^ap  "  ^ap  "  ''"'cx3  "  ^'oa  "^  STI  (e+n^)aVh) ,   we  can  reduce  the 
equations  (9.1)  to  two  equations  for  two  unknown  functions  * 

10 
and  03. 

For  notational  cowenience  we  define  the  invariant  dif- 

4 
ferential  operators   O   and   A 

o\u.v)  =  u|°v|^  -  u|°v|^  =  oNv.u) 

Under  the  assumption  that  the  shell  is  shallow,  that  the 
surface  traction  forces  can  be  derived  from  a  potential  function, 
and  the  body  forces  in  the  plane  of  the  shell  are  negligible, 
we  have  from  (9.1) 

(5a)   |[A(h'^Z^)  -  (l+v)C^(h"^,^)]  +  C^(f,'-o)  +  io^(a>,a)) 

-  |(|(S+  +  S-)  +  (l-v)|)  =  0'(a2(e+n2)/x2), 


For  an  engineer's  approach  to  this  problem,  together  with 
some  series  solution  to  the  resultant  eauations,  see 
Mansfield  [5]. 
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(s  +s~ )  I 

(5b) ^rkA(h^AD)  -  (l-v)0^(h^,a))]  +  0^(<:>,a>+f)  +  h  °^^  °^ 

.  SA(f4<D) 2^ 2_A(f  -  if^o.)  +  |-(S^  -  S-) 

where  Q  =  2h/x,   E  =  Young's  modulus,   v  =  Pols son's  ratio. 
CO  =  X  -f,   X  =  coordinate  of  the  deformed  midsurface,   f  =  the 
X-^  coordinate  of  the  undeformed  midsurface,   h  =  the  variation 
in  thickness  of  the  undeformed  shell,   <t>  =  the  Airy  stress 
function,   S.  =  the  surface  loads  in  the  U   coordinate  system 
per  unit  deformed  surface  area  on  the  deformed  shell,   S  =  a 
potential  function  such  that  S  ^  =  |-(S^  -  S^),   and   B,  =  body 
force  per  unit  deformed  volume  in  the  IP  coordinate  system. 

Note  that  by  our  assumptions  on  the  magnitude  of  the 
applied  forces,  the  forces  normal  to  the  midsurface  can  be 
of  order  e.   The  resultant  normal  force,  however,  must  be 
considerably  smaller.   In  other  words  we  have 

S^,  S^  =  £r(e), 

S3  -  S^  =  £f(en^). 
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APPENDIX 

General  Inequalities  for  Hilbert  Norms 
of  Functions  on  a  Slab 

These  results  are  independent  of  the  theory  of  elasticity. 
Greek  subscripts   a,^,...   range  over  the  values   1,2. 
Repeated  subscripts  or  superscripts  are  to  be  summed. 

Let   X   and   H  be  fixed  positive  quantities.   The 
function  w  =  w(U"'',U^,U^)   is  defined  and  of  class   C^   in 
the  region  MtU^U^'  <  X^,   |U^|  <  H. 

We  define  the  symbols   "g"   and   "o"   as  follows: 

A  =  [7(B)   with  B  >  0  means  that  there  exists  a  uni- 
versal constant   K  such  that   |A|  <  KB. 

A  =  ^(B)  +  c(C)   with   B  >  0  and   C  >  0  means  there 
exists  a  universal  function  K(6)   such  that  for  each  positive 
5  the  inequality 

|A|  <  K(5)B  +  5C 

is  satisfied. 

We  define: 

I'l^ii^  =  JjlwI^dU^dU^dU^  =  JlwI^dU, 
M  M 

w  =  sup  iw| , 

M 

3W 
w    - 


'^    oU^' 


AW  =  w^^^  +  w^22- 
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We   employ  an  avociliaiy  function 

$(u\u2)  =  (1  -  x-^\j%°-f. 

Clearly  for  U  =  (U-'-,u'^,U^)   in  M 

(Ao)     0  <  $  <  1,  ./r;^  <  2x-]/T7  v'^;^;^^  ^  ^^'^' 

1>      and  its  first  derivatives  vanish  for  U  U*^  =  x  , 
enabling  us  to  carry  out  integration  by  parts  with  respect 
to  U   without  producing  boundary  terms. 

Proofs  of  the  following  inequalities  are  based  on 
integration  by  parts. 


(Al)  il0  w  ^!!  =  &(x"^!!wll)  +  a(I!0awl!) 

PROOF 

,dU 


a,  3  M 


,ap  ,a  ,p  ,a  ,a  ,^\i 

M 

M 

+  °(J'(^^a*^a''^^3'''^3  +  l0Aw|2)dU) 
W      '       '       '       ' 

=  6"(x"^'>f )  +  a(  5^  "*,a'',3"^  +  "*AwH^). 


LSMivIA 


(A2)  !'#W  ^/  =  fTC^^Aw;;  +  \"^"w'i). 
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PROOF 

I   ll*",ae"'  =  I*'",ae",a3<^" 
a,p  M 

M 

a, 3  a, 3 

Using  Al  we  obtain  A2. 

LEMMA 

(A3)       "0w  „ll  =  ^((l+k)x"^'lwll)  +  o(k-^xr.0Awl'), 

where   k  is_  en  arbitrary  positive  number. 

PROOF 

2 


a  M 


=   -J(^   WW  ^^   +   20W  ^<p   ^w)dU 


r-/^2 

,  aa  >  ct   ,  a 

M 

=   0^((k^+l)x"^l!w:!^)    +  a(k"^X^:'^   wl'^   +  ^  i|0w      "^). 

a 


LEMMA 


(A^)  !'**', a^,p!!   =  ©"(wdl^AwJl   +    X"^!!w|!)). 


PROOF 


y  P"  «^  fll!^   =    [^^w     w     w  _w  „dU 
Z-  jOt   ,pi'  J         ,a   ,a   ,p    ,p 

a,p  M 


M 

=  gr(w2(ll0Awl|2   +   x-^llw!l2))    +0(7   l|*w     w  ,l|2), 
.         80  "'^ 


where  we  have  used  (AO)  and  (A2). 

More  complicated  to  derive,  are  estimates  involving  derivatives 

3  12 

with  respect  to  U  .   Let  h(U  ,U  )   be  defined  and  of  class 

C   in  the  region  U  U°^  <  >,  ,   and  such  that 
and  let 


H=   sup  h(U-'-,U^),   H*  =   inf  hCu-'-^U^),   w  ,  =  ^ -^. 


LEMMAS 


(A6)      P^a!^^\y\   =  &(X-^iiw|l)  +  a(x"^!Ih20w^^^l!), 

(A7)      !l0hw^^l!  =  9(liwl!)  +  adl^h^w^,^!!), 

(A8)      {Hhw^^^\\,      \\{^hw)^^J,      !!0(hw)^,^!|,   Il0hw^^^!!) 

=  0(x"^!;w!!  +  Ilh^w  ^jtl  +  H!l(/>Awj|)  +  a(x"0Awll) 

(A9)      1I0W  ,w  ,hl!  =  £r(w(ll0w  ,,hll  +  (H*)"^llwi|)). 

PROOFS 

In  the  proofs  of  all  these  inequalities  we  make  use  of 

an  extension  W  of  the  function  w  suggested  for  this  pur- 

1   2   "^j 
pose  by  L.  Nirenberg.   The  extension  W(U  ,U  ,U'^)   is  defined 

in  the  region  M-'-:U°U°  <  X^,   |  U^  |  <  2h  by 


W  =< 


3w(U-'-,U^,2H-U^)  -  2w(U-*-,U^,3H-2U-^)     for   H  <  U^  <  2H, 
w(u)  for   -H  <  U^  <  H, 

J'w(U-'-,U^,-2H-U^)  -  2w(U-'-,U^,-3H-2U^)   for   -2H  <  U^  <  -H. 


^ 


Clearly  W  and  W  -,      are  continuous  and  W  ^^  has  a  jump 
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■5 

discontinuity  at     U     =  ±H.      Moreover 


J 


2H 


-2H 


.S; 


— "^nF 


dU- 


■  »(C 


3^w 


(SU^ 


dU^j   for  k  =  0,1,2 


In  analogy  to  0  we  introduce  the  function 

3  3^ 


f(U^)  =   1  - 


^7" 


Then  following  the  proof  of  (A3), 


H 


2H 


2H 


r  h'^iw  ,|^du^  =  r^"  h^f^lw  ,|^du^  =  &((i+k^)r   Iwi^du^) 


-H 


,3 


-2H 


2H 


+  a(k"V;    f^lW  „l^dU^) 


J  o,r  '-,33 


-2H 

H  H 

=  C?((i+k2)r  Iwi^du^)  +  o(k-Vr  iw  ,,|2du^. 


-H 


Setting  k.  =  1  multiplying  by  h   and  integrating  over 

U-^,U^,   establishes  (A?).   Multiplying  by  h^^  jp      ,      putting 

,  a  ,  u 

2      2         1  12 

k  =  (>.  0  „0  „)"   and  integrating  over  U  ,U  ,   and  noting 
,  ct  ,  ex 

that  0  ^  =  C^(x'^)   gives  (A6). 

To  establish  (A8)  we  observe  that 

!!0hw^^^'l  =  er(i'(0fw)^^^hii  +  ii0^^fw^,hn 


+  "0f  ,w  h!!  +  11$  ^f  ,whll) 

=  ?r(|ih(0fw)  ,„!!  +  li*  w  ,h||  +  "0w  n  +  !!:/>  wH) 


>>■ 


a' 


82 


Now 


•J   V   tJ 

M 


2r    5    ../v,-! 


SU' 


•H(h"-^0wf) 


,a 


-1 


SU 


1      2      "5 
dU  dU  dU^ 


2„2r 


<    JJo^    H 


M' 


OU^SU'' 


(h      0Wf) 


12      3 
dU  dU  dU^ 


< 


2  2 

M 
r  f  r2hh  „H^— ^ — -(h"^$Wf  )-^(h"-'-^Wf  )dU-'-dU^dU^ 

Ml 


5U^  ^U 


BU' 


< 


< 


^Th^lh'-'-^Wf)        h  $fW  ,,+20f  ,W  ,+0f  ^,W  dU'^'dU^dU^ 
ojo  , eta  )P^  tP    >^        >^^ 

M^ 
2rj|[h(*Wf)    ^jrh  a(0Wf)    ^IdU^dU^dU^ 

jr;h2(h-i0Wf)^^^h0W^33fduidu2dU^ 
m1 


2/^-1 


lo,r2.„2 


2rrrh''(h"-'0Wf)    ^„irh^(0f  ^W  ^+0f  ^,W)]dU-'dU^dU 


M' 


+  l(c'!h(0fw)^3^l 


2   +  C-l|10h^^w^^f 


-^, 


-4, 


<  Cr(H^[!!'^>',aa''      +   '^        lU*^),all     "^   ^      !!*^li   +   x'^(  llwf +!|0wf )  ] 


+   "h^w  ,,!|2   +   x"^'l0hw  ,||2   +   X'^ll^w"^) 


,33 


,3' 


+  a(X(!lh(0w)^^^f  +   X^!!(*w),aaf  +   lK^^),af  ^   "^   >^"'^"*^ 


a 
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^iHiere  we  have  used  (AO),  (A5)  and  Schwarz's  inequality, 
and  C  is  a  constant.  Using  {kj>) ,  (A6)  and  (A7),  we  obtain 
(A8). 

Finally,  for  the  proof  of  (A9)  we  have 


j\2|w_,iW  =  0(j'"hV|w^,lW) 


-H     '^  ^  -2H 

-2H 


w2j   (fV|w^^3|2  +  h-2|w|2)du^) 


where  v/e  have  followed  the  reasoning  of  (A^l).   Multiplying  by 

2  12 

(0)    and  integrating  over   U  ,U   gives  (A9). 


SOBOLEV  S  INEQUALITY  FOR  A  SLAB.   For  U   in  M  we  have 

(AlO)  v^x"-'-l$(U)w(U)|  =  &(N!l0Aw"  +  Nx'^'iw"  +  ^"h^^w  ,,!'), 

X  N     '^^ 


where   N  >\A+(H>r. 


PROOF 


■5 
Assume  without  restriction  -H  <  U-^  <  0.   Let,  for  s  >  0, 


|V^|  <  H  if  s  <  Nx  <  -^  and   Hx"'''  <  1.   As  a  function  of 


V^   =  U^  +   s-^     where      '"^'^^   =  1     and      0  <   -^   <  H(Nx)'''".      Now 

H 

s,      V^(0)    =  U-^(O)      and 

IV^(Nx)    -   V^(0)|    <  H. 
But, 


Vv^  >/(v^-u")(v"-u^)   -^/^  >  s/ff  -   X 

>  Nx/l-c/e^    -    X    >   X(V^-H^X"^    -    !)• 
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Hence  if   n"^  >  4  +  H^x"^,   then  V^V^  >  X.   Denote  by  S   the 
lue  of  s   such  that  V^{s)V^{s)    =  >^.   Clearly  0  <  S  <  Nx. 


va 


Since  i),1>         vanish  for  s  >  S  we  can  write: 
>  ^ 

9  9  /  ^   H^  \' 

^^(U)w  (U)  =  (  r  s-^-^(U+SF)w(U+SP)ds  ) 
^  S    ^2       i  k    " 

=   r     ,   ,,(j>w)c^^^sds 


0  P^U  ?^ 


^^^ 


<  N^r 


0 


'  2 , 

s  ds 


3ince   F^'-"-!  =  -"-^  <  H'^(Nx)"^   ?^ 
0^(u)w2(u)  ^  ^^(NxrY^  |(?)w) 


are  bounded,  we  have 


aP> 


a,p 


+  h^Nx)"^^|(^w)^^^l2  +  h\Nx)"^*w^,^|2)s2ds 
a 


Now  integrating  over  all  unit  vectors   f  such  that 
0  <  p-^  <  H(Nx)~'^  produces: 


2ttH(Nx)"-^?)^(U)w^(u)  = 


a 


From  (Al),  (A2),  (A3),  (A5),  (A6),  (A?),  (A8),  follows  (AlO) 
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